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HIGHLIGHT  SUMMARY 

The  performance  of  fine  grained  soils  is  iiifluenced  by 
their  interaction  with  water.  The  permeability  and  the 
frost  heave  rate  are  among  the  properties  which  are  affected 
by  the  soil— water  interaction.  In  the  past,  various 
predictive  equations  based  on  volumetric  and  pore  size 
distribution  parameters  have  been  proposed  to  estimate 
permeability  and  frost  heave  rate.  However,  these  equations 
did  not  consider  the  effect  of  soil— water  interaction.  This 
study  investigates  this  effect  by  comparing  the  soil 
permeability  to  air  and  water. 

Another  concern  tc  engineers  is  the  ability  to  predict 
the  equilibrium  water  content  beneath  pavements.  A 
knowledge  of  the  so i 1  — moi sture  characteristic  curve  plays  a 
key  role  in  this  prediction.  Conventional  methods  to  obtain 
the  so i 1 —moi sture  characteristic  curve  involve  measurements 
of  soil  suction  and  are  time  consuming.  In  this  study  it 
was  shown  that  pore-size  distribution  data  can  be  used  to 
accurately  predict  the  practical  soi 1-moistui-e 
characteristic  curve  in  much  less  time.  Thus,  prediction  of 
in— place  equilibrium  v/ater  content  can  now  be  accomplished  on  a 
routine  basis. 


Natural  soils  and  the  artificial  mixtures  of  silt  and 
kaolin  were  used  in  this  study.  The  natural  soils  were 
compacted  by  impact  compaction  whereas  the  silt-kaolin 
mixtures  were  compacted  by  kneciding  compaction.  The  natural 
soils  had  been  obtained  by  co-operation  of  the  researchers 
at  the  University  of  Illinois. 

The  mercury  intrusion  technique  was  used  to  obtain  the 
pore  size  distribution  of  the  soils.  The  freeze  drying 
method  was  used  to  dehydrate  the  samples  p^rior  to  pore  size 
measurements . 

A  simple  failing  head  air  permeability  apparatus  was 
used  to  m.easure  the  air  permeability  of  the  samples.  Water 
permeability  values  were  taken  from  the  results  of  Garcia- 
Bengochea( 1978) .  So i 1 -moi sture  characteristic  curves 
obtained  by  conventional  method  were  provided  by  Janssen  and 
Dempsey( 1980) . 

The  pore  size  distribution  of  the  soils  tested  were 
predominantly  bi-modal.  However,  about  one— third  of  the 
samples  tested  exhibited  single  modal  characteristic. 
Varying  the  compaction  variables  produced  changes  in  the 
position  and  magnitude  of  the  pore  modes,  and  these 
variations  were  reflected  in  the  permeability. 


Two  permeability  models  which  relate  pore  size 
distribution  parameters  to  permeability  were  examined. 
These   models,   with   modifications,   were    successful    in 
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predicting  the  air  permeability,  but  they  failed  to  predict 
the  v/ater  permeability.  This  is  a  clear  indication  that 
minerology  has  an  influence  on  permeability.  The  clay 
content  and  composition  appear  to  be  the  controlling 
factors.  A  correlation  between  the  air  to  water 
permeability  ratio  and  soil  Activity  was  accomplished  to 
account  for  the  minerolcgic  effect. 


It  was  shown  that  the  minero 1 ogica 1  effects  on  clays  oV 

low   plasticity   are   very  low   and   hence,   the   predictive 

equations   developed   for  frost   heave   rate   in  such  soils 
require  no  modifications. 


INTRODUCTION 

The  performance  of  fine-grained  soils  is  influenced  by 
their  interaction  with  water.  Depending  on  the  "boundary 
conditions  and  purpose  of  the  project,  it  may  be  water 
flowing  that  is  the  concern.  At  other  times  it  may  be  how 
fast  the  water  will  freeze.  An  additional  concern  is  to 
what  equilibrium  magnitude  will  the  water  content  approach, 
and  this  will  influence  the  strength  and  deformation 
characteristics  in  service.  These  questions  have  been 
addressed  by  soil  engineers  for  a  long  time. 

The  above  matters  appear  to  be  different,  but,  in  fact, 
they  are  all  related  to  the  sizes  of  the  voids  within  the 
soils  and  their  continuity.  The  determination  of  these 
sizes  of  voids  has  been  difficult  in  the  past.  As  a 
consequence,  empirical  predictive  equations  have  been 
developed  based  on  gross  volumetric  parameters  or  grain  size 
distribution,  and  they  have  had  limited  success  for  sand. 

These  procedures,  however,  failed  for  the  fine-grained 
soils.  In  recent  years,  efforts  have  been  made  to  directly 
measure  the  pore  size  distribution  of  such  soils  and  relate 
these  to  the  soil  properties.  These  studies  produced 
suitable   prediction   procedures   for  soils  of  low  activity; 


therefore,  the  effect  of  minerology  was  not  considered. 

One  of  the  purposes  of  this  project  is  the  examination 
of  the  effect  of  soil  minerology  upon  the  prediction  of 
water  permeability  for  the  soils.  Minerology  has  an  effect 
by  "bonding"  water  to  the  sides  of  the  void  spaces  (by 
electro-chemical  attraction  forces)  to  change  the 
"effective"  pore  dimension  through  which  water  is  moving. 
In  the  present  study  an  attempt  has  been  made  to  shed  some 
light  on  the  effect  of  minerology  and  fabric  on  permeability 
and  freezing  rate  of  fine  grained  soils  by  comparing  their 
permeability  to  air  and  water. 

A  second  objective  of  the  study  is  to  develop  a  method 
to  predict  the  equilibrium  moisture  content  beneath 
pavements.  It  is  known  that  such  equilibria  are  reached 
under  certain  conditions.  Existing  methods  utilize  the 
soil-moisture  characteristic  curve  and  the  depth  to  water 
table  to  predict  the  moisture  content. 


Conventional  methods  to  obtain  the  soil-moisture 
characteristic  curve  involve  measurements  of  soil  suction 
and  are  time  consuming.  This  study  focuses  attention  on  the 
use  of  the  mercury  intrusion  technique  to  obtain  the  soil- 
moisture  characteristic  curve  by  relating  suction  to  pore 
size  distribution  characteristics.  This  procedure  is  much 
quicker,  and,  once  again,  recognizes  the  dependence  of  this 
behavior  upon  the  pore  sizes  which  control  it. 


The  mercury  intrusion  technique  was  used  to  obtain  the 
pore  size  distribution  of  the  soils.  The  freeze  drying 
method  was  used  to  dry  the  samples  before  intruding  with 
mercury . 

A  simple  falling  head  air  permeability  apparatus  was 
used  to  measure  the  air  permeability  of  the  soil  samples. 
Water  permeability  values  were  taken  from  the  results  of 
Garcia-Bengochea( 1978) .  Soil-moisture  characteristic  curves 
obtained  by  conventional  methods  were  provided  by  Janssen 
and  Dempsey(l980) . 


1  -  LITERATURE  REVIEW  AND  GENERAL  DISCUSSION 

1-1   Permeability  and  Fabric  of  Finegrained  Soils 

Permeability  is  one  of  the  important  properties  of 
soils  of  interest  to  geotechnical  engineers.  It  influences, 
among  other  things,  rate  of  settlement,  seepage,  dewatering, 
and  drainage. 

Permeability  of  compacted  soils  is  usually  measured  in 
the  laboratory  using  the  laboratory  compacted  samples.  When 
these  values  are  used  in  the  analysis  of  a  field  problem 
enormous  error  may  result.  The  laain  reason  for  this  is  the 
difficulty  in  replicating  the  field  soil  fabric. 

Pore  size  distribution  measurements  have  been  used  to 
predict  the  permeability  and  compare  the  fabric  of 
laboratory  and  field  compacted  samples.  One  of  the  main 
objectives  of  this  study  is  to  explore  the  effect  of 
minerology  on  permeability  by  comparing  the  air  and  water 
permeability.  and  to  develop  predictive  equations  for 
permeability  using  pore  size  distribution  measurements. 

This  section  discusses  the  importance  of  minerology   on 
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the   flow   of  water,  reviews  the  results  of  previous  studies 
concerned  with  the   permeability   and   fabric   of   compacted 


clays,  and  examines  two  permeability  prediction  models  which 
incorporate  pore  size  distribution  parameters. 

1-1.1   Effect  of  Minerology  on  Permeability 

Clay  minerals  are  very  small  colloidal  sized 
crystalline  partic 1 es ( diameter  less  than  1  micron),  derived 
primarily  from  chemical  weathering  of  certain  rock  forming 
minerals.  The  individual  crystals  look  like  tiny  plates  or 
flakes,  and  these  flakes  are  found  to  consist  of  many 
crystal  sheets  which  have  a  repeating  atomic  structure.  Two 
fundamental  crystal  sheets,  the  tetrahedral  or  silica,  and 
the  octahedral  or  alumina,  are  found  to  exist.  The 
particular  way  in  which  these  sheets  are  stacked,  together 
with  different  bonding  and  different  metallic  ions  in  the 
crystal  lattice,  constitute  the  different  clay  minerals.  It 
is  rather  common  that  the  cations  in  the  octahedral  sheets 
are  substituted  by  different  ions  of  approximately  the  same 
physical  size;  this  process  is  called  isomorphous 
substituion(Yong  and  Warkentin,  1975). 

Isomorphous  substitution  and  imperfections  in  the 
crystal  lattice  lead  to  negative  electric  charges  on  clay 
particles.  Cations  from  the  pore  water  are  attracted  to  the 
particles  to  maintain  electro-neutrality.  These  ions  are 
called  the  exchangeable  cations  and  their  number  is  the 
cation-exchange-capacity  or  the  amount  of  negative  charge 
per  unit  weight  or  per  unit  surface  area  of  the  clay. 
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Clay  particles  in  soils  are  always  hydrated,  i.e 
surrounded  by  layers  of  water  molecules  called  adsorbed 
water.  The  forces  holding  water  molecules  to  the  clay  arise 
from  hydrogen  bonding,  and  e 1 ectro— static  attraction.  The 
water  molecules  are  held  against  the  clay  mineral  surface  in 
many  layers  and  the  bonding  forces  become  weaker  with 
distance  from  the  surface.  The  density  and  viscosity  of 
this  water  close  to  the  clay  surface  is  higher  than  that  of 
Cree  water  (Low,  1961,  1979).  There  is  lack  of  agreement  as 
to  the  thickness  of  the  hydration  layer,  and  also  there  is 
no  sharp  division  between  water  of  hydration  and  free  water 
in  a  clay— water  system. 

In  a  dry  clay,  adsorbed  cations  are  tightly  held 
against  the  surface.  When  the  clay  is  placed  in  water, 
there  is  a  tendency  for  the  adsorbed  cations  to  diffuse  away 
from  the  clay  surface  in  order  to  equalize  the  concentration 
throughout.  Their  freedom  to  do  so,  however,  is  restricted 
by  the  negative  electric  field  originating  in  the  particle 
surfaces.  The  negative  surface  and  the  distributed  charge 
in  the  adjacent  phase  are  together  termed  the  diffuse  double 
layer.  The  interaction  of  diffuse  ion  layers  of  adjacent 
particles  is  given  as  an  explanation  for  swelling  of   clays. 

The  flow  of  water  through  clayey  soils  is  affected  by 
the  clay-water  interaction.  High  viscosity  of  water  near 
the  clay  surface  and  countere 1 ectro— osmosis  are  said  to 
affect  the  permeability.   The  other   factors   affecting   the 


flow  of  water  are  discussed  in  the  following  section.  If 
high  viscosity  exists  in  an  appreciable  volume  of  the  clay 
pore  space,  flow  rates  will  be  less  than  those  predicted 
from  Poisui lie's  law  which  assumes  that  the  pore  liquid 
viscosity  is  a  constant  throughout  and  equal  to  the 
viscosity  of  free    water. 

The  phenomenon  of  countere 1 ectro— osmosis  can  be 
explained  as  follows:  the  movement  of  water  past  the  clay 
surface  causes  a  downstream  transport  of  the  ionic  charges 
in  the  diffuse  part  of  the  double  layer;  this  displacement 
develops  a  streaming  potential  which  causes  e 1 ectro— osmosis 
in  a  direction  opposite  to  the  forward  flow;  the  net  result 
is  a  reduction  of  the  flow  rate  below  that  predicted  from 
Poisui lie's  1  aw. 

1-1.2   Permeability 

The  permeability  of  soils  has  been  of  interest  to 
researchers  for  a  long  time.  Darcy,  in  1856,  proposed  his 
well  known  law  for  flow  of  fluids  through  porous  media. 
Kozeny(1927)  suggested  an  equation  to  predict  the 
coefficient  of  permeability  using  the  physical 
characteristics  of  the  porous  medium.  Carman(1956)  modified 
Kozeny's  equation  and  since  then  it  has  been  called  Kozeny- 
Carman  equation.  However  this  equation  did  not  consider  the 
interaction  between  the  pore  fluid  and  the  porous  medium. 
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Several  researchers  have  attempted  to  measure  the 
permeability  of  clays  to  both  polar  and  non-polar  fluids. 
They  found  that  the  permeability  to  polar  fluids  was  very 
much  lower  than  to  non-polar  fluids.  The  applicability  of 
Darcy's  law  and  the  Kozeny— Carman  equation  has  also  been 
investigated,  and  discrepancies  were  found.  Various 
explanations  were  offered  for  these  discrepancies,  and  since 
these  are  of  importance  to  the  present  study  they  will  be 
discussed  in  detail  below. 

Macey(1942)  reported  rates  of  flow  for  benzene  100,000 
to  1,000,000  times  greater  than  that  for  water  in  his  clay. 
He  believed  that  the  most  important  cause  for  the  lower 
permeability  in  water  is  the  anomolous  viscocity  of  the 
water  near  the  clay  surface. 

Michaels  and  Lin(1954,  1955)  showed  experimentally  thsut 
the  permeability  of  kaolinite  decreased  markedly  as  the 
polarity  of  the  permeating  fluid  increased  and  that  the  most 
important  factor  controlling  the  permeability  of  kaolinite 
was  the  degree  of  dispersion  of  the  kaolinite,  in  the 
presence  of  polar  fluids.  They  also  concluded  that  electro- 
osmotic  counter  flow  might  be  of  some  importance  when 
aqueous  solutions  were  used  as  permeants  as  adsorbed  liquid 
surface  films  appeared  to  have  little  effect  on  the 
permeabi 1 i ty . 


Reeve(1953)  measured  air  and  water  permeability  for 
various  soils.  The  clay  content  of  the  soils  he  tested 
varied  from  2.3  to  38  percent,  and  the  exchangab 1 e— sodium- 
percentage  ranged  between  0  and  75.  He  found  that  the  air 
to  water  permeability  ratio  varied  with  exchangab 1 e— sodium- 
percentage  as  well  as  with  clay  content,  and  the  variations 
ranged  from  1  to  10,000.  He  attributed  these  variations  to 
structural  instability  caused  by  def 1 occulat ion  and 
dispersion.  Brooks  and  Reeve(1959)  also  obtained  similar 
results.  However,  these  authors  did  not  correct  their  air 
permeability  values  for  probable  gas  slippage  which  would 
have  resulted  in  higher  values  for  measured  permeability. 

Corey(1957)  found  that  the  permeability  ratio  for  loamy 
sand  was  about  2,  even  after  correcting  for  gas  slippage, 
and  suggested  that  the  ratio  might  be  much  higher  for  clayey 
soi 1 s . 


Mesri  and  01son(1971)  conducted  one  dimensional 
consolidation  tests  on  kaolinite,  illite,  and  smectite. 
They  used  polar  and  non-polar  fluids  as  pore  fluids,  and 
calculated  the  coefficients  of  permeability  using  Terzaghi's 
theory  of  consolidation.  Some  of  their  findings  are  as 
f o 1 1 ows : 

1)  The  coefficients  of  permeability  were  found  to  be 
largest  for  non-polar  fluids,  smaller  for  polar 
fluids  of  low  dielectric  constant,  and  lowest  for 
water. 
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2)  At  the  same  void  ratio  the  kaolinite  was  shown  to 
be  200,000  times  more  pervious  than  smectite  to 
water;  illite  was  200  times  more  pervious  than 
smect  i te . 

3)  The  coefficient  of  permeability  to  non-polar  fluid 
was  the  same  for  all  three  clays  at  a  given  void 
ratio. 

They  were  of  the  opinion  that  the  effects  of  adsorbed 
water  was  inadequate  to  cause  these  differences,  and 
explained  the  cause  of  differences  as  follows:  the  clay 
particles  form  stable  aggregates  and  the  main  fluid  flow 
takes  place  in  the  large  voids  between  the  aggregates;  in  the 
presence  of  polar  fluids  the  aggregates  break  down  and  there 
results  a  dispersive  structure  with  low  permeability;  the 
degree  of  dispersion  depends  on  the  polarity  of  the  pore 
fluid  and  the  physico-chemical  properties  of  the  soil 
mineral s . 

The  Kozeny-Carman  equation  has  been  found  to  be 
unsuitable  for  the  prediction  of  the  permeability  of  fine 
grained  soils  by  Lambe(1954),  Michaels  and  Lin(l954),  and 
01sen(1962).  Lambe(1954)  found  that  the  changes  in 
permeability  with  water  content  could  not  be  accounted  for 
by  changes  in  dry  density  or  void  ratio,  and  he  concluded 
that  they  must  have  been  caused  by  the  changes  in  the 
structure  of  the  soil  particles,  i.e.,  the  soil  fabric. 
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01sen(1962)  measured  hydraulic  flow  rates,  electrical 
conductivities,  and  streaming  potentials  for  flow  through 
kaolinite,  illite,  and  Boston  blue  clay.  He  found  discrep- 
ancies between  measured  flow  rates  and  those  calculated  from 
Kozeny-Carman  equation.  He  concluded  that  electro-osmotic 
counter  flow,  high  viscosities,  and  tortuous  flow  paths 
failed  to  account  for  the  permeability  characteristics  he 
measured,  and  suggested  that  unequal  pore  sizes  was  the 
cause  for  the  discrepancies. 

Schmid(1957)  suggested  that  an  adsorbed  layer  can  be 
considered  as  a  stationary  boundary  layer,  and  derived  an 
equation  for  permeability  using  effective  porosity. 

Non-Darcy  flow  has  been  observed  by  many  researchers. 
Von  Engelhardt  and  Tunn(1955)  measured  flow  rates  through 
clay-bearing  sandstones  and  found  that  the  flow  rates 
increased  more  rapidly  with  hydraulic  gradient  than 
predicted  by  Darcy's  equation.  They  interpreted  the 
behavior  as  resulting  from  the  water  behaving  as  a  non- 
Newtonian  liquid  whose  viscosity  depends  on  the  shearing 
force,  and  hence,  on  the  gradient.  Kemper(1960)  postulated 
that,  such  effects  results  from  e 1 ectro— kinetic  coupling  or 
electro-osmotic  counter  flow. 

Hansbo(1960)  and  Martin(1962)  suggested  that  this  type 
of  behavior  could  have  resulted  from  particle  movements 
leading   to   reversible   void    plugging    and    unplugging. 
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Low(1961)  argued  that  that  the  existence  of  a  quasi- 
crystalline  adsorbed  water  structure  could  explain  non- 
1 ineari ty . 

Miller  and  Low(1963)  observed  a  thershold  gradient,  and 
suggested  that  the  presence  of  quasi— crystall ine  water  and 
unequal  pore  sizes  were  the  reasons  for  this  behavior. 
01sen(1965)  suggested  that  undetected  experimental  errors 
could  have  resulted  in  the  apparant  threshold  gradient  in 
the  published  results  and  favored  a  particle  migration 
concept.  This  was  also  supported  by  Mitchell  and 
Younger(1967) . 

Miller  et  al.(1969)  reported  that  no  threshold  gradient 
was  found  for  soils  with  montmori 1 1 oni te  content  of  upto  50 
percent.  Russel  et  al.(1971)  also  found  no  evidence  for 
non-linear  flow  in  non— swelling  soil,  but  they  reported  non- 
linear flow  in  a  soil  having  small  content  of  swelling  clay. 

The  flow  of  gases  through  porous  media  has  been 
investigated  by  many  researchers  including  Muskat ( 1946) , 
Sullivan  and  Hertal ( 1940) ,  and  Kl inkenberg( 1941 ) .  Darcy's 
law  has  been  found  to  be  applicable  at  low  pressure 
gradient,  and  equations  were  derived  for  compressible  flow. 
It  was  noted  that,  with  samples  of  lower  permeability,  the 
values  obtained  with  gases  are  generally  higher  than  those 
obtained  using  liquids. 
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Adzumi (1937)  suggested  that  slip  occured  with  gases. 
Ddlcli  (1956)  explained  tlris  as  follnws:  "Wlien  tlu-  mean  froe 
path  of  gas  becomes  appreciable  with  respect  to  the  size  of 
the  channels  through  which  it  is  flowing,  the  viscous-flow 
assumption  of  zero  gas  velocity  at  the  channel  wall  no 
longer  holds  and  gas  begins  to  slip  past  the  channel  walls. 
The  slip  contributes  an  additional  component  to  the  velocity 
beyond  that  of  viscous  flow  and  this  component  becomes 
appreciable  at  low  pressures  and/or  when  the  channels  are 
extremely  small.  K 1 inkenbergC 1941 )  suggested  that  the  true 
permeability  can  be  obtained  by  measuring  the  apparent 
permeability  as  a  function  of  reciprocal  mean  pressure  and 
determining  the  intercept  of  the  plot. 

Several  researchers  including  Ladd(1960)  and  Langfelder 
et  al.(1968)  measured  air  permeability  of  partially 
saturated  kneading  compacted  clays,  and  reported  that  near 
optimum  water  content  air  permeability  was  essentially  zero. 
Blight(1971)  conducted  air  permeability  tests  on  soils  from 
tailings  deposits,  which  were  subsequently  compacted  and  air 
dried  in  the  laboratory,  and  found  non— linear  flow 
characteristics.  Serracin ( 1980)  measured  air  permeability 
of  compacted  grundite  and  found  that  Darcy's  law  was  obeyed. 

Summarizing,  recent  investigations  suggest  that  no 
threshold  gradient  exists  in  non— swelling  soils,  and  non- 
linearity,  if  found,  could  have  resulted  from  particle 
migration.   The  large  difference  in   permeability   to   polar 
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and   non-polar   fluids  is  due  to  the  break  down  of  aggregate 
structure  in  the  presence  of  polar  fluids. 

1-1.3   Fabric 

The  terms  fabric  and  structure  of  soil,  as  used  in  this 
study,  were  defined  by  Yong  and  Sheeran  (1973)  as  follows: 
"The  structure  of  the  soil  has  been  defined  as  that 
property  of  soil  which  provides  for  its  integrity.  An 
important  component  of  structure  is  its  fabric,  i.e., 
the  physical  arrangement  of  soil  particles  including 
particle  spacing  and  pore  size  distribution.  When  soil 
fabric  is  considered  in  conjunction  with  bonding  forces 
and  particle  interaction  mechanism  developed,  the 
structure  of  soil  is  obtained". 

The  fabric  can  be  divided  broadly  into  two  categories: 
macrofabric,  and  microfabric.  Macrofabric  can  be  seen  with 
the  naked  eyes  or  a  hand  lens,  and  it  includes  features  such 
as  stratification,  fissures,  voids,  and  large  scale 
inhoraogenei t ies .  Microfabric  requires  at  least  an  optical 
microscope  to  be  seen  (Mitchell,  1976).  Of  primary  concern 
to  this  study  of  clayey  soils  is  the  microfabric,  and  it  is 
this  category  of  fabric  that  is  discussed  in  more  detail 
be  1 ow. 

Several  conceptual  model  of  the  microfabric  of  natural 
soils  have  been  proposed  in  the  past,  the  earliest  of   which 
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were  conceived  by  Terzaghi ( 1925) ,  Go  1 dschimdt ( 1926) ,  and 
Casagrande( 1932) .  Terzaghi 's  model  assumes  that  the 
individual  clay  minerals  stick  to  each  other  at  the  points 
of  contact  with  forces  sufficiently  strong  for  the  creation 
of  a  "honeycomb"  arrangement.  Go  1 dschimdt ( 1926)  suggested 
that  the  clay  particles  in  highly  sensitive  clays  were 
arranged  in  an  unstable  card  house  arrangement. 
CasagrandeC 1932)  postulated  that  "honeycomb"  arrangements 
existed  in  sensitive  marine  soils  and  that  these  were 
basically  similar  to  Terzaghi 's  concept;  in  addition,  he 
also  suggested  that  the  particle  arrangements  located  in  the 
smallest  gaps  between  adjacent  silt  particles  were  highly 
compressed  whereas  those  located  in  the  larger  spaces 
between  silt  particles  had  undergone  relatively  little 
compression.  Since  then,  there  have  been  many  other  models 
suggested  by  various  researchers.  Collins  and  McGown(1974) 
summarized  the  various  models  postulated  to  that  time,  and 
suggested  a  few  modifications  to  some  of  the  earlier 
concepts. 

Lambe  presented  an  explanation  of  the  microstructure  of 
inorganic  soil  in  1953,  and  subsequently  suggested  a  model 
of  the  structure  of  compacted  clay(1958).  He  used  colloidal 
and  crystal  chemistry  to  postulate  his  model,  and  individual 
particles  were  his  predominant  units.  The  Lambe  concept 
suggested  that  on  the  dry  side  of  optimum  water  content,  the 
moisture   deficiency   causes   the   shrinkage   of   the  double 
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layer,  and  allows  the  particles  to  approach  sufficiently 
close  for  edge  to  face  attractions  to  develop;  this  results 
in  a  flocculated  structure.  On  the  wet  side  of  optimum 
water  content  the  fully  developed  double  layers  prevent 
close  approach  and  a  face  to  face  arrangement  dominates; 
this  is  termed  a  dispersed  structure.  Seed  and  Chan(1959) 
related  this  concept  to  observed  engineering  behavior  of 
clays  compacted  dry  and  wet  of  optimum. 

Since  then  various  researchers  who  could  not  explain 
their  experimental  results  using  Lambe's  model  suggested 
modifications  to  his  model.  01sen(1962)  postulated  that 
clay  particles  may  occur  in  clusters,  and  these  clusters 
make  up  the  soil  structure.  It  is  interesting  to  note  that 
Carman,  in  1959,  presented  a  similar  concept  for  fine 
grained  soils.  Mesri  and  01son(l971)  also  suggested  a 
similar  model  for  clays.  However,  recent  studies  of  soil 
fabric  using  transmission  and  electron  microscopy  gave 
strong  evidence  for  the  need  for  modifications  to  existing 
models.  Sloane  and  Kel 1(1966)  found  no  evidence  for  edge- 
to-face  structure  for  kaolinite  samples  compacted  on  the  dry 
side  of  optimum,  but,  rather,  found  an  essentially  random 
arrangement  of  particles  or  domains.  Smart(1969)  defined 
the  term  domain  and  gave  a  classification  of  domains  and 
sub-domain  structures. 

Harden  and  Sides(1970)  observed  that  there  has  been  a 
marked  difference  in  structure  at  low  magnifications  between 
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samples  compacted  dry  of  optimum  and  those  compacted  wet  of 
optimum,  but  little  difference  was  found  at  high 
magnifications.  Harden  and  Sides  hypothesized  that  at  low 
compaction  moisture  contents,  the  macropeds  have  a  high 
strength  and  are  able  to  resist  the  compaction  pressures 
without  much  distortion.  Thus,  the  voids  between  these 
macropeds  are  of  the  same  macroscale  and  are  filled  with 
air.  As  the  compaction  moisture  content  is  increased,  the 
macropeds  become  wetter  and  consequently  weaker.  Thus, 
during  compaction,  the  macropeds  are  more  easily  distorted 
and  the  size  of  the  interpedal  air  voids  is  reduced. 
Further  increase  in  the  compaction  moisture  content  causes 
the  macropeds  to  become  even  weaker  and  during  compaction 
they  are  distorted  to  fill  the  macro  pores  which  tend  to 
disappear.  Thus,  it  is  impossible  to  distinguish  the 
individual  macropores  and  at  low  magnification  the  soil 
appears  homogeneous.  Independent  studies  by  Hodek(1972) 
suggested  a  similar  model.  Yong  and  Sheeran ( 1973)  provided 
a  scheme  of  classification  which  identifies  the  various 
observable  fabric  units  with  reference  to  the  mode  of  fabric 
viewing . 

Although  most  of  the  above  analysis  were  strictly 
qualitative  a  few  researchers  attempted  to  quantify  the 
fabric.  Mi tche 1 1 ( 1956)  used  optical  microscopy  to  quantify 
the  orientation  of  clay  particles  by  using  a  measure  of 
light  transmission  through  the  fabric.   The   measurement   of 
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pore  size  distribution  of  compacted  clays  remained  a 
neglected  subject  for  a  long  time.  A  number  of  researchers 
from  Purdue  University  attempted  to  use  pore  size 
distribution  as  a  descriptor  of  fabric  of   compacted   clays. 

Diamond( 1970 ,  1971)  used  scanning  electron  microscopy, 
x-ray  diffraction,  and  pore  size  distribution  methods  to 
study  the  fabric  of  impact  compacted  clays.  His  electron 
micrographs  of  clay  samples  appeared  to  be  similar  to  that 
of  Barden  and  Sides(1970).  X-ray  diffraction  measurements 
showed  no  evidence  of  increased  particle  orientation  with 
increasing  molding  water  content  as  suggested  by 
Lambe(1958).  Pore  size  distribution  measurements  indicated 
that  content  of  the  fine  pores  (inter  particle  spaces,  0.02— 
O.OSfim)  was  not  affected  by  the  compaction  water  content. 
He  also  suggested  that  the  size  of  a  domain  is  about  Bfim. 

Sridharan  et  al .  (1971)  performed  pore  size  distribu- 
tion measurements  of  statically  compacted  clays.  They  found 
that  volume  changes  due  to  compaction  at  a  given  water 
content  was  at  the  expense  of  larger  pores.  They  further 
stated  that  the  different  clays  appear  to  have 
characteristically  different  overall  features  in  their  pore 
size  curves  and  these  variations  are  not  simply  correlated 
with  differences  in  the  grain  size  distributions. 

Ahmed(197l)  found  no  significant  difference  in  the 
shapes   of   pore   size   distribution   curves   between   clays 
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prepared  by  impact  and  kneading  compaction  to  a  common  water 
content  and  density.  Electron  micrographs  of  compacted  clay 
samples  obtained  by  him  appeared  to  be  similar  to  that  of 
Diamond(1971) . 

Bhasin(1975)  observed  that  increasing  the  compactive 
effort  on  the  dry  side  of  optimum  for  a  given  soil  decreases 
the  total  porosity  and  diminishes  the  fraction  of  larger 
pores.  Increasing  the  compactive  effort  on  wet  side  of 
optimum  had  little  effect  on  either  the  total  porosity  or 
the  distribution  of  pore  sizes.  Garcia— Bengochea( 1978)  also 
stated  that  the  changes  in  soil  fabric  was  generally  at  the 
expense  of  the  larger  pores. 

The  concept  of  a  deformable  aggregate  model  is  strongly 
supported  by  the  above  described  reults. 

White(1980)  studied  the  fabric  of  laboratory  and  field 
compacted  soils.  He  concluded  that  fabric  was  significantly 
different  for  both  soils.  The  differences  were  found  to 
exist  at  all  water  contents  except  those  near  the  optimum 
water  content.  The  differences  near  the  optimum  water 
content  may  have  been  due  to  the  inability  to  accurately 
define  that  water  content,  especially  for  the  field 
relationships.  Not  only  were  laboratory  and  field  fabrics 
different,  but  also  different  field  equipments  create 
significantly  different  fabric;  however,  impact  and  kneading 
compaction   in  laboratory  produced  fabrics  not  significantly 
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different. 

1-1.4   Permeability  Models 

Attempts  have  been  made  in  the  past  to  correlate 
permeability  with  certain  measurable  properties  of  the  soils 
such  as  grain  size  distribution,  porosity,  internal  surface 
area,  etc.  Ha2en(1911)  proposed  an  empirical  equation  based 
on  grain  size  distribution  to  predict  permeability  of  sands. 
Various  other  empirical  equations  were  also  proposed.  Since 
these  correlations  did  not  consider  the  pore  geometry  they 
had  only  limited  success.  Efforts  to  refine  this  approach 
by  introducing  empirical  grain  shape  and  packing  parameters 
have  not  won  general  acceptance.  Seme  researchers  attempted 
to  correlate  permeability  with  porosity  or  void  ratio 
(Jacob, 1946,  Tay 1  or , 1948) .  However,  due  to  the  strong 
dependence  of  flow  rate  upon  width,  continuity,  shape,  and 
tortuosity  of  the  conducting  channels,  this  kind  of  approach 
will  not  yield  satisfactory  resul ts(Hi 1 1 e 1 ,  1980). 

Kozeny-Carrnan  theory,  which  is  based  on  the  concept  of 
hydraulic  radius,  enjoyed  popularity  for  a  long  time. 
However.  its  applicability  to  fine  grained  soils  was 
questioned  by  various  researchers  as  discussed  in  Section 
1-1.2. 

Numerous  theoretical  models  have  been  introduced  to 
represent   the   porous   media  by  a  set  of  relationships  that 
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are  amenable  to  mathematical  treatment.   Such  models  include 
capillaric,  parallel,  serial,  and  branching  models. 

Puree  1 1 ( 1949)  used  a  modified  capillary  model  to 
calculate  the  permeability  of  sandstone  from  pore  size 
distribution  measurements.  Chi  Ids  and  Co  1 1 is-George(1950) 
and  Marshal  1 ( 1958)  introduced  probability  into  the  capillary 
model  theory  to  improve  the  model. 

Scheidegger ( 1974)  comprehensively  discussed  the  complex 
nature  of  porous  media,  the  limitations  of  conventional 
models  for  predicting  permeability,  and  presented  several 
models. 

Garcia-Bengochea( 1978)  used  the  capillary  models  by 
Purcel 1 (1948)  and  Marshal  1 ( 1958) .  and  the  hydraulic  radius 
model,  with  few  modifications,  to  predict  the  permeability 
of  artificial  silty  clays.  Juang(1981)  proposed  a  new  model 
which  he  used  to  predict  the  permeability  of  sands  and  sand- 
clay  mixtures. 

The  main  objective  of  this  present  study  is  to  find  the 
effect  of  minerology  on  permeability.  Thus,  it  was  decided 
to  use  the  simple  capillary  models  used  by  Garcia- 
Bengochea( 1978) ;  this  also  faciliated  the  use  of  some  of  his 
data  for  this  present  study.  The  hydraulic  radius  model  was 
not  used  because  of  the  following  reasons: 

1)  Hydraulic  radius  model  assumes  uniform  pores(Carman , 
1959).    The   fine   grained  soils  used  in  this  study 
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are  shown  to  have  non-uniform  pores. 
2)  Garcia— Bengochea  found  that  capillary  models  give   a 
better  correlations  than  hydraulic  radius  models. 

1-1.4.1  Capillary  model.  The  derivation  of  this  model  is 
given  in  detail  by  Garcia— Bengochea( 1978) ,  and  only  some 
important  points  will  be  repeated  below. 

The  porous  medium  is  modeled  as  a  set  of  parallel 
cylindrical  capillaries  of  varying  pore  diameters  extending 
from  one  end  of  the  porous  medium  to  the  other.  Each 
capillary  is  assumed  to  have  a  constant  radius  throughout 
its  length.  The  frequency  of  occurence  of  each  capillary 
size  is  determined  from  the  pore  size  distribution  curve. 

The  permeability  is  given  by  the  following  equation, 

k  =  Tn_„/"  ^2f(5)d5  (1-1) 

32m 

Where   k  =  permeability 

T  =  unit  weight  of  fluid 

II   =  absolute  viscosity  of  the  fluid 

n  =  porosity 

6    =  diameter  of  the  pore 

f(£)  =  porosity  density  function 

„/"  r(5)do    =  1  (1-2) 

Replacing  the  permeability,  k,  by  the  physical 
permeability,    K,   and   using   the   relationship   k  =  KT/^ . 
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equation ( 1-1 )  becomes 

K  =  C^E(d^)T\  (1-3) 

where   C^   is  a  shape  factor  which  is  equal  to  1/32  for 
cylindrical  pores.   Also, 

E(d2)  =  „J~  6^r(S)dS  (1-4) 

There   are   several  shortcomings  to  these  equations,  as 
pointed  out  by  Scheidegger ( 1974) : 

1)  all  pores  are  assumed  to  go  in  one  direction. 

2)  tortuosity  is  not  considered. 

3)  the  shape  factor  is  considered  constant  for  all  pore 
sizes . 

4)  the   value   of   E(d2)   is  extremely  sensitive  to  the 
frequency  of  the  larger  pores. 

5)  the  model  is  an  over  simplification   of   the   porous 
medium. 

1-1.4.2  Modified  capillary  mode  1 (Marshal  1 , 1958) .  Childs  et 
al.(l950)  and  Marshal  1 ( 1958)  employed  another  form  of 
capillary  model  using  the  pore  size  distribution  to 
calculate  permeability.  They  considered  two  unit  cross- 
sections  of  a  porous  medium  placed  together  such  that  the 
pores  from  one  surface  are  randomly  connected  to  those  of 
the  other  surface.  Each  of  the  cross— sections  contain 
cylindrical  pores  of  varying  diameters.  The  pore  size 
distributions  for  the  cross-sections  are  assumed  to  be 
identical . 

It  can  be  shown  that, 
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^    -  Tn^  oJ"oJ"  52f(6i)f(6j)d6id6.  (1-5) 

32// 

where  k,T,n,M,f(5)  are  as  defined  above. 

5  =  smaller  of  the  pore  sizes  S^    and  6j. 

Rewriting  equation  (1-5)  in  terms  of  pore  diameter  and 
physical  permeability: 

K   =   C^a^n^  (1-6) 

where 

3^       =    oj°°o/"    52f(5i)f(6j)d6id6j.  (1-7) 

In  summary,  the  physical  permeability  can  be  described 
by: 

K  =  C3  PSP  (1-8) 

where    PSP  is  a  pore  size  parameter,  and 

Cg  is  a  shape  factor. 
Table  1  presents  a  summary  of  the  models  disscussed. 

1—2   Soi 1— Moisture  Properties 

Subgrade  moisture  is  an  important  variable  in  the 
design  and  performance  of  pavements.  The  moisture  content 
of  soils  has  been  related  to  California  Bearing  Ratio  (which 
is  a  measure  of  shear  strength),  the  dynamic  or  resilient 
modulus  (which  affects  the  deflection),  and  the  unsaturated 
hydraulic  conductivity  (which  controls  the  flow  of  water 
during  frost  heave) .  Even  though  the  moisture  content  of  a 
subgrade  is  known  at  compaction,  the  field  moisture  content 
can  and  will  be  different  after  some  time   in   service.    It 
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Table  1     SUMMARY  OF  PERMEABILITY  MODELS  USED 

1)  VARIABLE  DIAMETER  CAPILLARY  MODEL 

Equation  K  =  C^E(d2)n 

Pore  Size  Distribution 

Parameter  E(d2)n  =  ^f'^d^V  (  S)  dS 

2)  MARSHAL  MODEL 

Equation  K  =  C^H^n^ 

Pore  Size  Distribution 

Parameter  3^  =  ^  J~^  J~S2  f^  (  5  j  )  f  (  5  j  )  d6  i  d5  j 

where  S  is  smaller  of  the 
pore  sizes  5i  and  5 j 


K  =  Physical  Permeability,  L^ 

n  =  porosity 

C3  =  shape  factor 

5i  =  pore  diameter,  L 

f(5)  =  porosity  density  function,  „f°°    f(5)d(5  =  1 
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would  be  very  useful  to  have  a  rational  method  to  predict 
the  long— term  moisture  distribution  beneath  pavements.  In 
the  following  section  the  concept  of  matric  potential  or 
soil  suction,  and  the  soi 1— moisture  characteristic  curve  is 
explained;  in  addition  an  existing  rational  method  to 
predict  the  moisture  content  is  discussed. 

1-2.1   Soil-Moisture  Potential 

Soil  water  contains  energy  in  both  potential  and 
kinetic  forms,  just  as  any  other  bodies  in  nature.  As  the 
water  moves  very  slowly  in  soil,  the  kinetic  energy 
component  is  usually  ignored  in  the  study  of  the  soil— water 
system.  The  potential  energy  of  the  soil  water  can  be 
subdivided  into  gravitational  potential,  pressure  potential, 
osmotic  potential,  and  gas  potential.  The  gravitational  and 
pressure  potential  are  the  primary  force  fields  in  the  soil- 
water  system.  The  osmotic  potential,  which  results  from  the 
presence  of  solute  in  the  system,  and  gas  potential,  which 
is  dependent  upon  external  or  internal  gas  pressure,  are 
considered  to  have  minor  influence  on  the  total  potential, 
and  are  normally  ignored. 


The  total  potential  can  be  written  as  follows: 

H  =  Z  +  h  (1-9) 

where         Z  =  gravitation  potential  or  position  head 

h  =  pressure  head  or  matric  potential 

H  =  total  head 
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When  there  is  no  moisture  movement,  total  head  is  equal 
to  zero . 

Therefore , 

h   =  -Z  (1-10) 

The  matric  potential  results  from  capillary  and 
adsorptive  forces  developed  in  the  soil  matrix.  It  can  be 
visualized  as  a  suction  exerted  by  the  soil  to  pull  water  up 
from  the  water  table.  This  suction  is  dependent  on  the  soil 
pore  geometry  or  soil  matrix;  thus  the  term  matric  potential 
is  used(Hil lei  ,  1980) . 

The  concept  of  matric  potential  can  be  easily 
understood  by  an  anology  made  between  the  pores  in  the  soil 
and  capillary  tubes.  Water  in  a  narrow  tube  will  rise  above 
the  free  water  level  until  an  equilibrium  condition  is 
reached.  The  water  in  the  tube  forms  a  meniscus  at  the  air- 
water  interface;  the  dimensions  of  the  meniscus  are  related 
to  the  size  of  the  tube.  The  smaller  the  radius  of 
curvature  of  the  meniscus,  the  higher  the  water  will  rise  in 
a  capillary  tube.  A  similar  meniscus  can  form  in  the  water 
between  soil  particles.  The  height  of  capillary  rise  can  be 
described  by  the  equation. 


h  =  4  T  Cose 
Td 


(1-11) 


where,  h  =  height  of  capillary  rise 
T  =  surface  tension  of  water 
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6   =  contact  angle  between  the  wall  of   the   tube 

and  water. 
T  =  unit  weight  of  water 
d  =  diameter  of  the  tube 

1-2.2   Soi 1— Moisture  Characteristic  Curve 

A  graph  showing  the  relationship  between  moisture 
content  and  suction  is  called  the  soil-moisture 
characteristic  curve.  Since  experience  shows  that  soil 
suction  can  take  on  a  large  range  of  values,  it  is  common  to 
plot  the  logarithm  of  suction  in  cm  of  water  against 
moisture  content  (Figure  1).  Several  soil  scientists 
including  Keen(1931)  and  Scho'f ie  1  d(  1935)  used  this  curve  to 
study  some  aspects  of  the  soil  structure.  The  nature  of  the 
capillary  pressures  and  the  role  they  play  in  petroleum 
reservoir  behavior  have  been  discussed  by  Levrett ( 1941 )  and 
others . 

Several  researchers  including  Croney(1952)  of  the 
British  Road  Research  Laboratory,  have  used  suction-moisture 
content  relationships  to  predict  moisture  distribution 
beneath  road  and  airfield  foundations.  as  well  as  road 
performance.  Green  and  Corey(1971)  used  this  relationship 
to  predict  unsaturated  hydraulic  conductivity.  Elzeftawy 
and  Dempsey( 1977)  used  soil  moisture  characteristic  data 
along  with  measurement  of  saturated  hydraulic  conductivities 
to   predict  unsaturated  hydraulic  conductivities  for  highway 
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Figure   I     Soil  Moisture  Cliaracteristic  Curve 
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soils.  The  overall  goal  has  been  to  formulate  a  theory 
which  could  reliably  predict  the  moisture  movement  beneath 
pavements,  because  moisture  movements  affect  the  strength 
parameters  considerably. 

Croney  et  al.(1951)  have  described  the  methods  used  to 
determine  the  soil  moisture  characteristic  curve.  These 
methods  consist  of  the  tensiometer,  the  direct  suction 
method,  the  pressure  plate  method,  and  the  centrifuge 
method.  Usually  no  single  method  can  be  used  to  cover  the 
entire  moisture  tension  range,  so  several  methods  may  be 
used  in  actual  practice.  Puree  1 1 ( 1949)  used  mercury 
porosimeter  measurements  to  predict  the  suction  curve. 

The  soi 1 -moisture  characteristic  curve  can  be  obtained 
by  two  methods: 

1)  In  desorption,  by  taking  an  initially  saturated 
sample  and  applying  increasing  suction  to  gradually 
dry  the  soi 1 ,  and 

2)  In  sorption,  by  gradually  wetting  an  initially  dry 
samp  1 e . 

These  two  methods  usually  give  two  different  curves; 
this  is  the  well  known  hysteresis  effect.  According  to 
Hi  1 1  el (1980) ,  the  hysteresis  effect  is  attributed  to: 

1)  the  geometric  non— uni f ormi ty  of  the  individual 
pore,  called  the  "ink  bottle"  effect. 

2)  the   difference   in   contact  angle  of  the  advancing 
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and  receding  meniscus. 

3)  entrapped  air 

4)  swelling,  shrinking  or  aging  phenomena. 

The  amount  of  water  retained  at  low  values  of  matric 
suet ion(between  0  and  1  bar)  depends  primarily  upon  the 
capillary  effect,  and  hence  is  strongly  affected  by  the  soil 
structure.  On  the  other  hand,  water  retention  in  the  higher 
suction  range  is  due  increasingly  to  adsorption  and  is  thus 
influenced  less  by  the  structure  and  more  by  the  texture  and 
specific  surface  of  the  soil  material. 

1-2.3   Prediction  of  Soil-Moisture  Characteristic  Curve 

Purcel 1 (1949)  suggested  a  method  to  get  the  soil- 
moisture  characteristic  curve  from  mercury  intrusion 
measurements.  His  method  is  as  follows:  In  the  conventional 
method  to  get  the  soil-moisture  characteristic  curve  water 
is  used  as  the  pore  fluid.  Since  it  wets  the  soil 
surface(i.e.  contact  angle  is  less  than  90°),  surface  forces 
are  induced  and  these  forces  draw  water  into  the  soil  pores. 
However,  removal  of  water  from  pores  will  require  external 
pressure.  On  the  other  hand  mercury,  a  non  wetting  fluid, 
requires  external  pressure  to  intrude  the  pores.  Since  the 
Washburn  equation  applies  to  both  cases,  if  the  measurements 
with  respect  to  mercury  are  known,  then  the  capillary 
pressures  for  water  can  be  predicted,  provided  the  surface 
tension  and  the  contact  angle  for  both  fluids  are  known. 
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Washburn's  equation  is 


P„  =  4  T„  Cos0, 
d 


(1-12) 


where   Tj 
al  SO  , 


surface  tension  of  mercury, 
pressure  required  for  mercury, 
contact  angle  between  mercury  and  soil 


P...  =  4  T...  Cose. 


(1-13) 


where   T^^  =  surface  tension  of  water 

P^  =  pressure  required  for  water. 

6^   =  contact  angle  between  water  and  soil 

From  equation  (1-12)  and  (1-13), 

T   Cos0 


(1-14) 


The  soil-moisture  characteristic  curve  predicted  from 
mercury  intrusion  measurements  is  equivalent  to  the 
desorption  curve.  At  each  intrusion  pressure  to  obtain  the 
corresponding  moisture  content,  the  intruded  volume  of 
mercury  should  be  subtracted  from  the  total  volume  of  pores. 

The  primary  advantage  of  this  method  is  that  it  allows 
one  to  obtain  the  soi 1 -moisture  characteristic  curve  over  a 
wide  range  of  pressures  within  an  hour  whereas  the 
conventional  method  takes  weeks  and  also  requires  the  use  of 
different  equipment. 
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1-2.4   Prediction  of  In-Situ  Equilibrium  Water  Content 

Croney  and  Co  1 eman ( 1953)  used  the  matric  potential 
(which  they  called  the  pore  water  pressure)  and  soil- 
moisture  characteristic  curve  to  predict  equilibrium 
moisture  contents  below  road  surfaces,  where  the  ground 
water  table  is  at  shallow  depths(i.e.  about  20  feet  in  clay 
soils). 

The  principles  on  which  their  method  is  based  are  as 
f o 1 1 ows : 

1)  The  pore  water  pressure(matr ic  potential)  tends  to 
an  equilibrium  value  dependent  only  on  the  position 
of  the  ground  water  table. 

2)  There  is  a  relationship  between  pore  water  pressure 
in  the  soil  at  a  given  location  and  the  suction  of 
this  soil,  where  suction  is  the  pore  water  pressure 
of  a  sample. 

3)  There  is  a  relationship  between  suction  and  water 
content(soil  moisture  characteristic  curve). 

Equilibrium  conditions  are  as  follows: 

a)  The  subgrade  soil  is  at  a  constant,  uniform 
temperature  above  0°  C,  and 

b)  The  subgrade  soil  cannot  imbibe  or  lose  moisture 
through  the  pavement  or  adjacent  soil  masses. 

At  equilibrium  the  pore  water  pressure  or  matric 
potential      can  be   expressed  by   the   following   equation 
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u    =    -Z  (1-15) 

where      u      is    the    pore    water   pressure    expressed    in    cm   of 

water   and   Z   is    the   height   above   the   ground  water   table.        Suction    is 
taken   to   be   the   value   of   pore  water  pressure   of   an  undisturbed   sample 
which,    owing   to   its   removal,    has  been  relieved   from  all   external   stresses. 
The   suction  of   a   soil   can  be   represented  by   the   following   equation, 

s   =  u  -  ap  (1-16) 

where,    p    is    the    overburden    pressure. 

a.      is      the    compressibility    factor    and    it    varies 
between    0    and    1. 

It  has  been  demonstrated  that  oc  is  given  by  the 
derivative  of  the  shrinkage  curve  of  the  soil(Croney  et 
al. ,1958).  If  there  is  no  equipment  available  to  obtain  a 
shrinkage    curve    the    following    relationship    may    be    used, 

a.   =    0.027PI    -    0.12  (1-17) 

where,  PI  is  the  plasticity  index  of  the  soil. 

The  moisture  content  can  be  directly  read  off  from  the 
soil-moisture  characteristic  curve  using  the  value  of 
suction  calculated. 
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1—3  Pore  Size  Distribution 

The  measurement  of  pore  size  distribution  allows  the 
engineer  to  obtain  a  quantitative  measure  of  the  fabric  of 
soils.  It  has  frequently  been  noted  that  the  grain  size 
criteria  are  insufficient  to  predict  the  behavior  of  soils 
under  varied  conditions.  Since  it  is  the  pore  sizes  which 
affect  the  permeability  and  change  during  compression, 
sample  disturbance,  and  swelling,  the  study  of  pore  size 
distribution  is  an  important  tool  to  the  engineer. 

1-3.1  Applications 

Pore  size  distribution  studies  have  been  of  interest  to 
researchers  in  various  fields  including  chemical 
engineering,  petroleum  engineering,  soil  science,  and 
geotechnical  engineering.  It  has  been  used  to  predict  the 
permeability  by  Puree  1 1 ( 1949) ,  Childs  et  al.(1950), 
Marshal  1 (1958) ,  Mi  1 1 ington  and  Quirk(1959),  Green  and 
Corey(1971),  Garc ia-Bengochea( 1978) ,  Juang(1981),  and 
others . 

Si  1 veria(1965)  combined  pore  size  distribution 
calculations  with  probabilistic  modeling  to  study  the 
problem  of  fine  particles  passing  through  protective  sand 
filters.  Sridharan(1968)  used  pore  size  distribution 
studies  to  calculate  the  negative  pore  pressures  in  partly 
saturated   clays  in  an  effort  to  better  understand  the  shear 
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strength  of  partly  saturated  soils   by   examining   the   soil 
structure . 

Ahmad(1971)  related  pore  size  distribution  to  the 
strength  of  clayey  soils.  Reed(1977)  correlated  frost  heave 
rate  of  compacted  soils  with  pore  size  distribution.  The 
application  of  pore  size  distribution  in  many  branches  of 
engineering  has  been  comprehensively  discussed  by 
Ahmad(1971)  and  Garc ia-Bengochea( 1978) . 

Several  methods  are  available  to  measure  the  pore  size 
distribution.  The  gas  adsorption  technique,  capillary  rise 
test,  pressure  plate  suction  test,  and  mercury  intrusion  are 
among  the  commonly  used  methods.  Bhasin(1975)  presented  a 
critical  review  of  the  various  techniques  available  and 
concluded  that  only  the  mercury  intrusion  method  is  capable 
of  measuring  the  entire  range  of  pore  sizes  of  fine  grained 
sols.  The  mercury  intrusion  method  is  used  in  this  study 
following  the  lead  of  previous  researchers  of  Purdue 
University.       , 

1—3.2   Theory,  Assumptions,  and  Limitations 

The  basis  of  the  mercury  intrusion  method,  as 
originally  conceived  by  Washburn( 1921 ) ,  and  subsequently 
applied  by  Ritter  and  Drake(l945),  rests  on  the  principle 
that  a  non-wetting  fluid  will  require  external  pressure  to 
intrude   a   porous  medium.   The  pressure  required  depends  on 
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the  contact  angle,  the  pore  shape,  and  the  surface  tension 
of  the  liquid.  For  cylindrical  pores  Washburn  suggested  the 
following  relation: 

(1-18) 

where,  P  =  required  pressure 

T   =   surface   tension   of   the  intruding 

fluid 
a   =    contact  angle  between  the   solid   and 

the  fluid 
d  =  diameter  of  the  pore 

Mercury  has  been  found  to  be  the  best  intruding  fluid 
because  of  its  low  vapour  pressure,  inertness,  and  non 
wetting  property  to  most  surf aces(Wins 1 ow  and  Diamond, 1970) . 
Diamond( 1970)  measured  the  contact  angle  for  mercury  on 
various  minerals  by  the  sessile  drop  method.  He  suggested 
that  contact  angle  on  kaolinite  or  illite  was  within  one 
degree  of  147  degrees,  and  on  montraori 1 1 oni te  was  within  one 
degree  of  139  degrees. 

With  the  assumption  of  constancy  of  the  surface  tension 
and  contact  angle,  Washburn's  equation  implies  that  as 
pressure  increases  mercury  wi 1 1  be  intruded  into 
progressively  smaller  pores.  For  each  pressure  increment, 
mercury  is  forced  into  the  accessible  soil  pores  of  a 
diameter  larger  than  or  equal  to  that  calculated  from 
equation  (1-18)  at   the   corresponding   pressure.    Intruded 
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volume  is  recorded  along  with  the  pressure  as  the  experiment 
progresses . 

The  most  important  limitation  of  this  method  is  that 
the  pore  diameter  measured  at  a  given  pressure  is  the 
diameter  of  the  channel  leading  into  the  pore,  which  could 
be  smaller  than  the  actual  size  of  the  pore(Ritter  and 
Drake,  1945).  As  the  result  of  this  "ink  bottle"  effect, 
the  pore  diameters  measured  are  refered  to  as  "limiting  pore 
diameters".  Cebeci  et  al.(1978)  attempted  to  evaluate  the 
arrangement  of  ink  bottle  pores  by  a  second  intrusion 
technique.  However,  obtaining  the  true  pore  size 
distribution  is  still  a  difficult  task.  This  problem  is  of 
minor  importance  in  the  case  of  f 1 ov/  through  porous  media, 
because  the  flow  channels  through  which  the  water  has  to 
flow  are  probably  equivalent  to  chose  intruded   by   mercury. 

Other   assumptions   or   factors   to   be   considered  for 
mercury  intrusion,  as  listed  by  Garcia-Bengochea( 1978) ,  are: 

1)  Only   the   pores  accessible  from  the  out  side  of  the 
sample  are  penetrated. 

2)  The  minimum  pore  size  is  dependent  on   the   pressure 
capacity  of  the  porosimeter. 

3)  The   soil   is  considered  to  be  incompressible  during 
pressur i zat ion  and  does  not  break  down. 

4)  The  contact  angle  and   surface   tension   values   are 
assumed  to  be  constant  during  pressuri zat ion . 

5)  The   mercury   will   not   react   with  the  soil  during 
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intrusion . 
6)  Cylindrical  pore  shapes   are   assumed   to   calculate 
pore  diameter. 

A  detailed  discussion  of  the  assumption  and  sources  of 
error  of  the  mercury  intrusion  technique  can  be  found  in 
Ritter  and  Drake(l945),  Rootare( 1968) ,  Wins  1 ow( 1969) , 
Diamond(1970) ,  Bhasin( 1975) ,  Garcia-Bengochea( 1978) ,  Kenney 
(1980).  Reed  et  al.(1980).  and  Juang(1981). 

1—4   Soil  Drying  Methods 

The  measurement  of  pore  size  distribution  of  a  soil  by 
mercury  intrusion  technique  requires  dry  samples.  The 
dehydration  process  may  cause  shrinkage,  expansion, 
cracking,  fabric  reorientation  or  even  minero 1 ogical 
changes.  Several  methods  have  been  developed  to  dehydrate 
the  samples.  The  amount  of  change  in  gross  volume  has  been 
the  most  common  criterion  used  to  monitor  the  fabric 
changes.  Minero 1 ogica 1  changes  can  be  detected  by 
differential  thermal  analysis.  X-ray  diffraction,  and 
scanning  electron  microscopy. 

•Air  or  oven  drying  is  the  simplest  method  of 
dehydration  and  was  used  by  Diamond( 1970)  and  Sridharan  et 
al.(1971)  for  a  compacted  clay.  The  authors  concluded  that 
oven  drying  method  did  not  affect  their  results  adversely. 
However,  they  noted  air  or  oven  drying  is  not   an   effective 
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drying  method  for  many  saturated  or  non— saturated  clayey 
soils,  because  of  the  large  shrinkage  that  usually  occurs 
with  drying. 

Critical  region  drying  is  the  most  complex  procedure  of 
dehydrating  soils,  and  was  discussed  by  Bhasin( 1975) .  The 
method  consists  of  elevating  the  temperature  and  pressure  of 
a  saturated  soil  specimen  in  a  controlled  manner  to  the 
critical  region  of  water,  at  which  the  gas  and  liquid  exist 
in  a  single  phase.  Once  the  critical  region  is  reached,  the 
pressure  and  temperature  of  the  sample  are  slowly  reduced  to 
atmospheric  in  a  way  that  water  in  the  sample  remains  in  the 
gas  phase.  Thus,  the  shrinkage  forces  caused  by  air— water 
menisci  are  eliminated.  However,  the  minerology  of  the  soil 
could  be  changed  by  the  high  temperatures  and  pressures  used 
to  reach  the  critical  region  of  water. 

The  fluid  replacement  method  is  another  drying 
technique  discussed  by  Tovey  and  Yan(1973).  This  method 
involves  replacing  the  pore  water  with  another  liquid  of  low 
surface  tension,  usually  an  organic  liquid.  This  liquid  is 
removed  by  air  or  oven  drying.  Since  the  shrinkage  forces 
are  small  because  of  the  low  surface  tension,  there  will  not 
be  substantial  change  in  the  soil  fabric.  However,  it  was 
found  that  swelling  frequently  accompanies  fluid 
replacement.  For  this  reason  and  because  of  the 
complication  of  the  process.  fluid  replacement  is  not 
considered  suitable  for  large  investigations. 
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Freeze  drying  is  the  method  used  in  this  study  as 
recommended  by  Zimmie  and  Alraal eh ( 1976) .  This  method 
consists  of  quickly  freezing  a  small  soil  specimen  to 
cryogenic  temperatures,  followed  by  vaccum  drying  the  frozen 
sample  to  remove  water  by  sublimation.  The  quick  freezing 
prevents  frost  heaving  by  minimizing  ice  crystal  nucleation 
and  water  migration.  Sublimation  eliminates  the  formation 
of  air-water  menisci,  and  hence  shrinkage  forces  which  alter 
the  soil  structure.  Zimmie  and  Alraal eh ( 1976)  pointed  out 
that  vaccum  drying  could  be  performed  at  room  temperature 
without  a  refrigerant,  because  sublimation  is  an  evaporative 
process  accompanied  by  cooling,  so,  although  the  specimen 
container  is  exposed  to  room  temperature  the  sample 
temperature  will  stay  below  freezing.  By  adopting  his 
recommendation  the  drying  technique  is  simplified  and 
accelerated. 

If  properly  conducted,  freeze  drying  of  soils  has  been 
shown  to  be  a  superior  method  for  removing  water  by  many 
researchers  including,  Diamond(1976) ,  Tovey  and  Yan(1973), 
and  Zimmie  and  A Ima 1 eh( 1976) .  The  volume  changes  are 
generally  less  than  5  percent.  For  this  reason  and  its 
relative  simplicity  freeze  drying  was  chosen  for  use  in  this 
study. 
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2  -  SOILS  STUDIED 


Two  types  of  soils  were  used  in  this  study:  an 
artificial  mixture  of  soil  materials,  and  number  of 
naturally  occuring  soils.  The  artificial  mixture  is 
composed  of  a  natural  silt  and  a  commercial  kaolin  clay. 
The  artificial  mixture  was  used  to  compare  with  the  watei: — 
permeability  values  obtained  by  Garc ia-Bengochea( 1978)  for 
the  same  materials.  These  soils  were  chosen  by  him  to 
represent  the  natural  silty  clays  which  are  frequently  used 
in  earthwork  construction. 

The  silt  was  taken  from  a  natural  loess  deposit  near 
the  Wabash  River  in  southwestern  Indiana.  The  soil  was 
excavated  from  the  second  bench  level  on  the  east  side  of 
U.S. Highway  41,  about  one  mile  south  of  Patoka,  Indiana. 
The  silt  was  air  dried  and  passed  through  a  No.  40  U.S 
standard  sieve  before  testing.  X-ray  diffraction  analysis 
showed  that  the  silt  was  composed  primarily  of  quartz  and 
dolomite  with  small  amounts  of  calcite  and  feldspar  also 
present.  The  silt  was  found  to  contain  1  to  3  percent  of 
clay  particles.  A  separate  analysis  on  the  clay  revealed 
the  presence  of  illite,  kaolinite,  and  raontmori 1 1 oni te  with 
a  possible  trace  of  vermicu 1 i te (Garc ia-Bengochea ,  1978). 
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The  kaolin  used  in  this  study  is  the  commercially 
available  Edgar  plastic  kaolin.  It  is  processed  to  remove 
material  coarser  than  40  fim  in  size,  dried  in  a  tunnel  dryer 
to  300°  to  450°F,  and  then  pulverized.  The  manufacturer 
states  that  99.5  percent  of  the  clay  mineral  composition  is 
kaolin  with  the  remainder  being  quartz  or  mica. 

Three  blends  of  silt  and  kaolin  were  tested: 
90?  silt  -  10^  kaolin 
70%  silt  -  30%  kaolin 
50%  silt  -  50%  kaolin 

Grain  size  distribution  curves  and  soil  classification 
properties  for  the  silt  and  for  each  of  these  mixtures  are 
shown  in  Figure  2  and  Table  2  respectively. 

The  natural  soils  used  in  this  study  were  taken  from 
various  places  in  Illinois,  and  provided  to  the  author  by 
Professor  W.J.Dempsey  of  the  University  of  Illinois.  A 
complete  description  of  these  soils  can  be  found  in  Janssen 
and  Dempsey ( 1980) .  The  relevant  properties  of  the  soils 
used  and  their  code  designations  are  shown  in  Table  3.  The 
soils  chosen  have  Activities  varying  from  0.13  to  2.2.  It 
was  hoped  that  the  variation  in  Activity  will  reflect  the 
various  degrees  of  soil— water  interaction  and  would  allow 
the  examination  of  a  much  larger  range  of  minerologic 
effects.  X-ray  diffraction  analyses  were  performed  on  most 
of  these  soils  and  the  results  are   given   in   Table  4.    As 
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Table  4   Results  of  X-ray  Diffraction  Analysis  on  Natural 
Soi 1 s  Tested 


Soil 


Minerals  Present 


ALVIN-B 

AVA-B 

CISNE-B 

HERRICK-A 

HOSMER-A 

HOSMER-B 

RICHVIEW-B 

SABLE-B 
STOY-A 


Illite,    Kaolinite,    Vermicu 1 i te ,    and 
random  mixed  layer  Vermiculite  and  Illite 

Kaolinite,  Illite,  Montmor i 1 1 oni te 

Kaolinite,  Montmori 1 1 oni te 


Kao linite.Illite, Montmori 1 1 onite , 
of  Vermiculite 


Trace 


Kaolinite,  Illite,  Vermiculite,  Random 
Mixed  Layer  Vermiculite  and  Illite 

Montmori 1 1 onite ,  small  amounts  of  Illite, 
Halloysite,  and  Kaolinite 

Kaolinite,  Illite,  Montmori 1 1 oni te , 
Random  Mixed  Layer  Vermiculite  and  Illite 

Kaolinite,  Strong  Montmori 1 1 oni te ,  Illite 

Kaolinite,  Illite,  Montmor i 1 1 oni te- 
disordered.  Mixed  Layer  Vermiculite  and 
II  lite 
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testing  progressed  it  was  found  that  the  results  of  X-ray 
diffraction  analysis  were  not  much  of  use  for  correlation 
purposes.  Consequently,  this  analysis  was  not  continued  for 
rest  of  the  soils. 
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3  -  APPARATUS  AND  EXPERIMENTAL  PROCEDURE 

3-1  Soil  Mixing  and  Curing 

The  mixing  of  silt  and  kaolin  was  achieved  using  a 
Patterson-Kel ley  liquid-solid  mechanical  blender.  The 
blender,  shown  in  Figure  3,  consists  of  a  rotating  V-shell 
inside  which  is  a  high  speed  rotating  dispersion  bar  from 
which  deionized  water  is  sprayed.  A  graduated  cylinder  is 
connected  to  the  dispersion  bar  through  a  stop  cock  which 
controls  the  amount  of  water  to  be  sprayed. 

To  use  the  blender,  2  Kg  of  the  air  dried  silt-clay 
combination  was  first  dry  mixed  for  about  10  minutes. 
Deionized  water  was  then  added  at  a  rate  of  50  to  75  ml  per 
minute.  After  the  desired  quantity  of  water  was  added,  the 
blender  was  stopped  and  soil  was  scrapped  off  the  inside 
walls  of  the  shell.  Following  scraping,  the  mixer  was  run 
for  about  7  minutes,  followed  again  by  scraping  and  another 
7  minutes  mixing  period.  The  soil  was  then  removed  from  the 
shell,  sealed  inside  a  polyethylene  bag,  and  cured  for  two 
days . 

When  adding  water  to  the  dry  mix,  experimentation 
indicated  that  adding  0.5  percent  less  water  than  calculated 
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Figure  3    Patterson-Keliey  Blender 
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gave  within  0.25  percent  of  the  desired  moisture  content. 
This  procedure  was  adopted  because  the  silt  and  kaolin  were 
at   some  smal 1  initial  water  content  in  the  air  dried  state. 

3-2  Compaction 

3-2.1  Apparatus 

An  electrically  driven,  semi-automatic  kneading 
compactor (Figure  4)  made  by  August  Manufacturing  Company  was 
used  to  compact  the  silt-kaolin  mixture.  The  sample  is 
compacted  by  a  triangular  shaped  foot  with  curvature  on 
sides.  The  pressure  applied  by  the  foot  is  controlled  by  a 
pneumatic-hydraulic  system  with  an  air  regulator  and  a 
pressure  gauge.  The  sample  holder  rotates  60  degrees 
between  each  tamp  in  order  to  achieve  uniform  compaction. 

The  pressure  applied  on  the  soil  by  the  foot  is  not 
directly  measured.  Gaudette( 1960)  and  Reed(1977)  have 
determined  calibration  curves  for  converting  gauge  pressures 
to  the  actual  foot  pressures  applied  to  the  soil.  However, 
Garcia-Bengochea( 1978)  suggested  that  these  curves  are  not 
unique  and  considerable  variations  in  the  foot  pressure  may 
occur  with  seasonal  variations  in  temperature  and  humidity. 
He  further  stated  that  if  the  foot  was  not  at  the  bottom  of 
its  stroke  before  tamping  started,  the  first  few  tamps  gave 
exceedingly    high    foot    pressures.     Bailey(1976)    and 
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Figure  4    Kneading  Compactor 
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Gaudette(1960)   give   further  details  about  the  operation  of 
the  kneading  compactor. 

Considerable  difficulties  were  encountered  in 
replicating  the  dry  densities  obtained  by  Garcia- 
BengocheaC 1978)  because  of  the  above  mentioned  problems. 

3-2.2  Procedure 

As  the  results  of  Garc ia-Bengochea( 1978)  were  used  in 
this  study,  the  procedure  adopted  by  him  to  compact  his 
samples  were  used.  The  samples  were  compacted  in  a  mold  4.0 
in.  in  diameter  and  1/30  cu. f t . (approx imate)  in  volume(ASTM 
D  698).  A  circular  filter  paper  of  diameter  4  in.  was 
placed  in  the  bottom  of  the  mold  to  prevent  the  soil  from 
adhering  to  the  base  plate.  The  soil  was  placed  and 
compacted  in  five  layers  of  approximately  equal  thickness. 
Thirty  tamps  of  the  compaction  foot  were  applied  to  each 
layer.  Before  adding  a  fresh  layer  of  soil  the  previously 
compacted  layer  was  scarified  in  order  to  prevent  separation 
at  the  contact  between  layers.  The  fresh  soil  was  tamped 
three  times  with  a  piston  before  compaction  began  in  order 
to  densify  the  soil  just  enough  to  prevent  the  mechanical 
foot  from  punching  through  the  layer. 

Once  the  compaction  procedure  was  completed,  the  collar 
was  removed  and  the  soil  was  trimmed  flush  with  the  top  of 
the  mold.   The  weight  of  the   mold   and   sample   was   noted. 
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Then  the  sample  was  extruded  from  the  mold,  and  cut  into 
small  cubes  for  freeze  drying-.  Trimmings  were  used  to 
obtain  the  water  content. 

Several  problems  as  noted  by  Garcia-Bengochea( 1978) 
were  encountered  during  compaction.  It  was  not  possible  to 
compact  the  soils  with  water  content  3  to  4  percent  greater 
than  the  optimum  water  content;  the  soil  tended  to  stick  to 
the  foot  during  its  upward  stroke,  preventing  effective 
compaction.  Also  when  the  moisture  content  was  slightly 
above  the  optimum  water  content,  the  soil  adjacent  to  the 
foot  heaved  and  cracked. 

The  natural  soils  obtained  from  Illinois  were  provided 
in  compacted  form.  The  soils  were  compacted  in  a 
commercially  available  Tempe  pressure  cell  (5.35  cm.  in  I.D 
and  3.0  cm.  in  height).  The  samples  were  prepared  by  mixing 
soil  and  water  to  give  the  optimum  water  content  for 
Standard  Proctor  compact ion(ASTM  D  698).  The  soil  was  then 
placed  into  a  Tempe  cell  and  compacted  by  20  blows  of  a  1.5 
lb  weight  falling  12  in.  The  sample  was  then  inverted  and 
compacted  an  additional  20  blows  on  the  other  side.  This 
procedure  gave  the  density  similar  to  that  specified  by  ASTM 
D  698.  A  detailed  description  of  the  procedure  can  be  found 
in  Janssen  and  Dempsey ( 1980) . 
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3-3  Freeze  Drying 

The  freeze  drying  procedure  as  suggested  by  Zimmie  and 
A 1  ma  1  eh ( 1976)  was  used  in  this  study. 

3—3.1  Apparatus 

The  apparatus  used  in  this  study  was  assembled  by 
Reed(1977)  and  the  following  description  is  taken  directly 
from  his  work.  Figures  5  and  6  show  the  basic  components  of 
the  freeze  drying  equipment  which  include 

1)  Wire  cage  to  hold  the  samples  during  freezing  and 
sub  1 iraat ion . 

2)  Large  Dewar  flask  containing  liquid  nitrogen  into 
which  the  sample  and  cage  are  dipped. 

3)  Dessicator  with  inside  wire  support  from  which 
sample  cage  is  suspended  during  sublimation. 

4)  Vaccura  pump  capable  of  evacuating  to  less  than  0.01 
mm  of  mercury. 

5)  Condenser,  which  is  connected  between  the  dessicator 
and  vaccum  pump  in  order  to  prevent  vapor  from  the 
soil  entering  the  pump.  Liquid  nitrogen  contained 
in  the  dewar  flask  was  used  to  cool   the   condenser. 

3-3.2  Procedure 

After  compaction  the  sample  was  extruded  from  the  mold 
using   an   hydraulic   extruder.    The   extruded   sample   was 
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Figure  5     Sample  Cage  Suspended  in  Desiccator 
(From  Reed,    I  977) 
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Figure  6    Freeze  Drying  Apparatus  with  Vacuum  Pump, 
Condenser  and  Desiccator  (From  Reed,   I  977) 
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dissected  into  several  pieces  using  a  wire  saw.   Small  cube  samples  8  mm 
on  a  side  were  trimmed  from  these  pieces  using  a  single  edge  razor  blade. 
The  trimmings  were  used  to  determine  the  water  content  of  the  sample. 
About  ten  small  specimens  were  cut  from  top,  bottom,  and  middle  portions 
of  the  sample.   The  trimmed  samples  were  placed  in  a  cage  and  lowered 
into  liquid  nitrogen.   They  remained  in  liquid  nitrogen  for  approximately 
10  min.  To  prevent  air  bubbles  from  adhering  to  the  sides  of  the 
specimens,  and  hence  insulating  the  samples  from  liquid  nitrogen,  the 
wire  cage  was  moved  up  and  down  during  its  stay  in  liquid  nitrogen. 

After  about  10  min.  the  vacuum  pump  was  started  and  the  wire  cage 
was  quickly  transferred  to  the  dessicator  for  vacuum  drying.   Vacuum 
was  maintained  at  approximately  0.1  mm  of  mercury  while  sublimation  took 
place.   Zimmie  and  Almaleh  (1976)  found  that  5  to  10  hr.  of  vacuum 
drying  would  be  sufficient  for  kaolinitic  soils,  and  9  to  10  hr.  would 
be  sufficient  for  montmorillonitic  soils.   Since  a  variety  of  soils  were 
used  in  this  study,  a  minimum  10  hr.  of  vacuum  drying  was  used  (usually 
the  vacuum  pump  was  run  overnight).   Once  sublimation  was  completed, 
the  samples  were  removed  from  the  wire  cage,  placed  in  labeled  bottles, 
and  stored  in  a  dessicator  containing  anhydrous  calcium  sulphate 
(commercial  name  "Drierite")  which  removed  any  remaining  moisture. 


59 


3-4  Pore  Size  Distribution 

The  apparatus  and  procedure  used  in  this  study  to 
determine  the  pore  size  distribution  is  essentially  the  same 
as  that  used  by  the  previous  Purdue  researchers  Ahmed(1977), 
Bhasin(1975) .  Reed(1977),  and  others.  Thus,  only  a  brief 
description  of  the  apparatus  and  procedure  will"  be  given 
be  1 ow. 

3-4. 1  Apparatus 

The  apparatus  used  in  the  measurement  of  pore  size 
distribution  is  composed  of  a  penetrometer;  a  filling 
device;  a  vacuum  pump;  a  McLeod  gauge;  a  mercury  manometer; 
a  control  board;  and  an  Aminco  porosimeter. 

The  penetrometer  (Figure  7)  consists  of  a  glass  tube 
attached  to  a  finely-calibrated  capillary  stem.  The  stem  is 
calibrated  in  0.002  ml  increments  and  has  a  capacity  of 
0.02  ml.  The  bulb  is  approximately  6.1  ml  in  capacity  and 
is  used  to  house  the  sample.  A  greased  stainless  steel  cap 
and  a  teflon  fastener  are  used  to  seal  the  bulb. 

The  filling  device  is  a  two  chambered  glass  tube  which 
holds  the  penetrometer  in  a  horizontal  position  during 
filling  and  low  pressure  intrusion.  The  filling  device  is 
connected  through  a  pressure  control  board  to  a  vacuum  pump, 
McLeod  gauge,  mercury  manometer,  and  a  bleeder  valve.  The 
McLeod  gauge  is  used  to   measure   pressures   below   1  mm   of 
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Figure  7    Penetrometer:  Disassembled  and 
Assembled  with  Sample 
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mercury  and  the  manometer  is  used  to  measure  pressures 
between  1  mm  of  mercury  and  atmospheric  pressure.  The 
arrangement  of  the  system  is  shown  in  Figure  8. 

An  Aminco  porosimeter (Figure  9)  is  used  to  apply 
pressures  above  atmospheric  pressure.  The  porosimeter  holds 
the  penetrometer  vertically  in  a  chamber,  and  applies 
pressures  up  to  15000  psi  with  an  electrically  driven 
hydraulic  pump  using  ethanol  as  the  pressurizing  fluid. 
Volume  changes  as  small  as  0.0001  ml  are  measured  using  an 
electro-mechanical  sensing   device  which  travels  in  the 

stem  of  the  penetrometer.  Two  Bourdon  pressure  gauges  of 
capacities  1000  psi  and  15000  psi  are  used  to  measure  the 
intruding  pressure. 

3-4.2  Procedure 

Following  freeze  drying,  the  samples  were  trimmed  with 
a  razor  blade  in  a  method  intended  to  produce  tensile 
failures  on  the  surface  to  keep  surface  pores  from  being 
smeared  shut.  After  trimming,  the  sample  was  weighed  to  the 
nearest  0.0001  g,  then  sealed  in  a  penetrometer  and  weighed 
again.  The  penetrometer  was  then  placed  in  the  filling 
device  and  sealed  around  the  stem  with  the  bulb  protruding. 
The  vacuum  pump  was  then  used  to  evacuate  the  system  to  less 
than  0.02  mm  of  mercury,  a  process  which  usually  took  30  to 
45  min.  • 
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Figure  8    Filling  Device,  McLeod  Gauge  and  Control  Board 
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Figure  9    Aminco  Porosimeter 
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After  evacuation  of  the  system  was  complete,  the 
penetrometer  was  filled  with  mercury.  This  was  accomplished 
by  rotating  the  filling  device  so  that  the  mercury  it 
contained  covered  the  opening  at  the  end  of  penetrometer 
stem.  The  vacuum  pump  was  then  turned  off,  air  was  bled 
into  the  system  to  a  pressure  of  approximately  20  ram  of 
mercury.  This  pressure  is  called  the  filling  pressure. 
Raising  the  pressure  forced  mercury  into  the  stem  and  bulb, 
and  surrounded  the  sample  to  be  intruded. 

After  the  mercury  completely  filled  the  penetrometer 
stem,  the  filling  device  was  rotated  again  back  to  its 
original  position.  This  broke  off  the  column  of  mercury 
near  the  end  of  the  stem  producing  an  initial  stem  reading 
which  was  recorded  along  with  the  filling  pressure.  This 
pressure  was  then  raised  in  increments  according  to  the 
equation  Pi/Pj-i  =  IC^,  where  c  was  0.18,  to  produce  a 
logarithmically  equal  class  width  as  suggested  by  Garcia- 
Bengochea(1978) . 

When  atmospheric  pressure  was  reached,  the  penetrometer 
was  removed  from  the  filling  device  and  weighed  in  order  to 
find  the  amount  of  mercury  in  the  stem.  A  final  stem 
reading  was  taken  with  the  axis  of  the  penetrometer 
horizontal;  it  was  then  placed  in  the  porosimeter.  An 
initial  reading  of  volume  intruded  was  taken  using  the 
electro-mechanical  probe  of  the  porosimeter.  These  steps 
allow  for  a  transition  from  low  pressure   to   high   pressure 
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readings.  The  pressure  was  increased  incrementally 
maintaining  intervals  as  for  the  low  pressure  readings  until 
the  pressure  was  just  below  the  capacity  of  15000  psi. 
Probe  readings  were  taken  after  each  increment  of  pressure 
was  added.  Only  a  short  time  (about  15  sec)  was  found  to  be 
necessary  to  complete  each  incremental  intrusion.  This 
process  (above  atmospheric)  was  known  as  high  pressure 
intrusion. 

3—5  Air  Permeability 

3—5. 1  Apparatus 

The  apparatus  (Figure  10)  used  in  this  study  to 
determine  air  permeability  was  developed  by  Serracin( 1980) 
and  modified  by  the  author  to  allow  for  the  use  of  a  drying 
mechanism.  The  schematic  diagram  of  the  apparatus  is  shown 
in  Figure  61,  Appendix  B.  "The  principle  of  the  apparatus  is 
that  the  flow  rate  is  measured  by  following  the  decay  of  a 
pressure  differential  imposed  across  the  sample.  A 
knowledge  of  the  decay  rate,  the  sample  dimensions,  and 
certain  apparaus  constants  permits  the  calculation  of  the 
permeability"  (Dolch,  1956). 

The  components  of  the  apparatus  are:  an  u-tube ,  one  end 
connected  to  a  glass  bulb  and  the  other  end  to  a   stopcock; 
a    sample    holder;    a    meter    stick;    and    a   plastic 
tube(commercial 1 y  known  as  dry  tube)  filled   with   a   drying 
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Figure   10    Air  Permeability  Apparatus 


67 

agent.    Mercury   or   water   is   used  as  the  manometer  fluid 
depending  on  the  permeability  of  the  sample. 

The  glass  bulb  at  one  end  of  the  u-tube  is  used  to 
support  the  sample  holder  during  the  testing.  It  is  tapered 
inside  to  facilitate  housing  the  sample  holder,  which  is 
also  tapered  to  fit  inside  the  bulb.  The  stopcock  is  used 
to  maintain  a  pressure  difference  before  starting  the  test. 
The  air  is  dried  by  allowing  it  to  pass  through  the  drying 
agent  before  it  enters  the  sample.  The  meter  stick  is  used 
to  measure  the  elevation  of  the  manometer  liquid  column 
during  the  testing. 

3-5.2  Procedure 

Following  freeze  drying,  the  dimensions  of  the  sample 
were  measured  using  a  vernier  caliper.  The  sample  holder 
was  placed  on  a  clean  glass  plate,  and  the  sample  was  placed 
in  position  in  the  sample  holder  and  sealed  in  place  by 
pouring  hot  paraffin.  Then  the  sample  holder  with  glass 
plate  was  placed  in  a  dessicator  until  the  paraffin 
hardened.  In  this  study,  it  was  found  that  applying  a  thin 
film  of  rubber  cement  to  the  sides  of  the  sample  and  the 
sample  holder  before  pouring  paraffin  gave  a  better  seal  and 
grip.  Also  to  avoid  excessive  shrinkage,  hot  paraffin  was 
poured  when  its  temperature  was  about  3°C  above  its  melting 
point . 
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In  the  meantime  a  pressure  difference  in  the  manometer  was  es- 
tablished by  applying  a  pressure  to  the  stopcock  side  of  the  u-tube  and 
maintained  by  closing  the  stopcock.   The  sample  holder  was  then  removed 
from  the  dessicator  and  a  thin  film  of  vacuum  grease  was  applied  to  the 
outside  before  placing  it  in  position  in  the  glass  bulb.   The  dry  tube 
was  connected  to  the  top  of  the  sample  holder  via  a  cork  and  a  tube. 

The  stopcock  was  then  opened  and  the  rate  of  movement  of  the 
manometer  level  was  recorded  using  a  stop  watch.   It  was  found  that  the 
time  required  for  each  test  varied  from  one  minute  to  one  hour  depending 
on  the  permeability.   The  equation  used  to  calculate  the  permeability 
and  its  derivation  are  given  in  Appendix  B. 
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4  -  RESULTS  AND  DISCUSSION  OF  RESULTS 

4—1  Compaction 

Three  levels  of  compaction  effort  were  used  for  the 
silt-kaolin  mixture.  The  medium  effort  was  chosen  by  the 
previous  researchers  to  obtain  an  optimum  dry  density  for 
90%  silt-  10%  kaolin  mixture  similar  to  that  obtained  by 
AASHTO  T99-70  (ASTM  D  698)  procedure( i . e .  Standard  Proctor) 
for  the  same  soil.  The  high  compaction  effort  was  used  to 
approximate  that  of  the  modified  Proctor  compaction 
procedure(ASTM  D  1557).  A  low  compaction  effort  was  also 
used. 

In  this  study , samples  of  the  silt-kaolin  mixture  were 
prepared  at  water  contents  selected  to  give  a  spread  in  the 
permeability  values.  It  was  not  possible  to  replicate  the 
dry  density  for  some*  samples  as  obtained  by  the  previous 
researchers  because  of  the  difficulties  discussed  in  Section 
3-2. 

The  code  designations  used  to  identify  the  silt-kaolin 
mixtures  are  given  in  Table  5.  The  code  gives  the  soil 
type,  the  compaction  effort,  and  water  content  relative  to 
its  optimum  water  content. 
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Table       5  Sample    Code    Designation    of      Artificial  Mixtures 

of    Silt   and   Kaolin 


Soil    Type 

S9  -  90^  Silt  and  10%  Kaolin 
S7  -  70%  Silt  and  30%  Kaolin 
S5  -  50%  Silt  and  50%  Kaolin 

Compaction  Effort  (Gauge  Pressure) 

L  -  Low  Compaction  Effort;  4.0  psi. 

M  -  Medium  Compaction  Effort;  8.5  psi. 

H  -  High  Compaction  Effort;  40  psi. 

Moisture  Condition 

0  -  Near  Optimum  Water  Content 
D  -  Dry  of  Optimum  Water  Content 
W  -  Wet  of  Optimum  Water  Content 

Examp 1 e 

S5MD:    50%    silt    and    50%    kaolin    compacted   at    8.5    psi 
gauge    pressure     and   dry  of   optimum  water   content. 
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All  natural  soils  used  in  this  study  obtained  from 
Illinois  were  compacted  at  their  respective  optimum  water 
content.  The  compaction  procedure  discussed  in  Section  3— 
2.2  was  used  to  approximate  the  maximum  dry  density  obtained 
by  the  procedure  of  ASTM  D  698  for  the  same  soil. 

Table  6  gives  the  sample  codes  and  the  compaction 
parameters  of  all  the  soils  tested.  Figures  11,  12,  and  13 
present  the  compaction  curves  for  silt— kaolin  mixtures  and 
the  location  of  each  of  the  sample  points.  The  compaction 
curves  were  obtained  by  Garcia-Bengochea( 1978) . 

4-2  Freeze  Drying 

The  freeze  drying  process  usea  in  this  study  was  simple 
and  suitable  for  drying  a  large  number  of  samples  at  one 
time.  Several  initial  measurements  indicated  that  volume 
change  during  freeze  drying  was  within  four  percent  as 
suggested  by  Garcia-Bengochea( 1978) ;  hence  no  additional 
checks  were  made  beyond  the  initial  ones. 

4-3  Permeabi 1 ity 

4-3.1  Air  Permeability 

The  air  permeability  was  measured  by  a  falling  head 
permeability  apparatus.  A  pressure  gradient  as  large  as 
4000   kN/m^   was   used.    The   calculation   procedures    are 


Table  6  Compaction  Parameters  of  Soils  Tested 


72 


SAMPLE 
CODE 

MOISTURE 

CONTENT 

(PERCENT) 

DRY  DENSITY 
(lb/ft3) 

S9HW 

13.5 

114.8 

S9MW 

16.8 

107.6 

S7H0 

12.97 

118.7 

S7M0 

14.3 

113.9 

S7L0 

18.5 

109.3 

S5M0 

17.7 

110.7 

S5MD 

18.5 

113.  1 

S5LD 

18.01 

109.0 

S5L0. 

19.6 

106.  1 

ALVIN-B 

12.8 

117.4 

AVA-B 

19.6 

106.7 

BRYCE-A 

24.4 

87.2 

CATLIN-B 

20.8 

85.7 

CISNE-B 

21.6 

91.8 

HERRICK-A 

22.4 

94.0 

HOSMER-A 

20.  1 

96.5 

HOSMER-B 

20.3 

96.6 

RICHVIEW-B 

19.6 

101.4 

SABLE-B 

21.6 

95.4 

STOY-A 

17.9 

98.4 
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Figure   I  I     Compaction  Curves  for  90%  Silt  and   10%  Kaolin 
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Figure    I  2    Compaction  Curves  for  70%  Silt  and  30%  Kaolin 
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Figure   I  3    Compaction  Curves  for  50%  Silt  and  50%  Kaolin 
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discussed  in  Appendix  B.  The  temperature  was  not  controlled 
during  the  testing;  it  was  monitored,  and  appropriate  values 
of  viscosity  and  unit  weight  of  air  at  that  temperature  were 
used  in  the  calculations. 

The  permeability  value  assigned  to  each  soil  is  the 
average  value  of  permeabilities  obtained  from  three  to  four 
samples.  For  each  sample  the  test  was  repeated  two  to  three 
times  and  the  average  value  was  taken  as  the  permeability  of 
the  sample.  When  the  permeability  test  was  repeated  for  the 
same  sample  the  values  of  permeability  obtained  were  in 
close  agreement  with  each  other,  whereas  among  samples  of 
the  same  soil  a  difference  of  one  order  of  magnitude  was 
found  in  some  cases. 

Figures  14  through  23  show  some  of  the  typical  plots  of 
flow  rate  versus  pressure  gradient.  It  is  clear  that 
Darcy's  law  is  obeyed  for  gases.  However,  a  small  intercept 
on  the  gradient  axis  was  obtained  in  many  cases.  Two 
possible  explanations  can  be  offered.  One  is  that  during 
testing  the  manometer  level  varied  very  rapidly  for  soils  of 
large  permeability  and,  hence,  the  recorded  readings  may  be 
slightly  in  error.  This  undei — estimation  of  the  rate  of 
flow  would  create  a  small  intercept  in  the  gradient  axis. 
The  other  possible  error  could  have  resulted  from  the  change 
in  the  contact  angle  of  the  manometer  fluid  on  glass  due  to 
atmospheric  contamination  of  the  glass  as  suggested  by 
01sen(1965).   It  was  found   that   for   a   0.635  cm   diameter 
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Figure  I  7   Air  Flow  Velocity  Vs   Pressure  Gradient  for  Stoy-A 
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Figure  I  9  Air  Flow  Velocity  Vs   Pressure  Gradient  for  S5L0 
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tube,  the  possible  error  in  the  gradient  assuming  contact 
angles  of  36°  and  76°  for  receding  and  advancing  meniscus 
respectively  was  about  3.21  KN/M^ .  The  combined  effects  of 
these  two  factors  could  have  yielded  a  small  intercept  on 
the  gradient  axis. 

It  was  also  found  that  drying  of  air  had  no  significant 
effect  on  the  permeability  of  the  soils.  In  order  to  check 
this,  some  samples  were  first  tested  with  dried  air.  Then 
the  samples  were  left  in  the  atmosphere  for  one  day  and 
tested  again.  As  can  be  seen  from  Figures  20  through  23  no 
significant  difference  in  permeability  was  found.  It  is 
also  possible  that  the  method  used  to  dry  the  air  was 
ineffective . 

The  physical  permeability  values  obtained  were  not 
corrected  for  probable  gas  slippage.  As  discussed  in 
Section  1-1.2,  gas  slippage  tends  to  give  higher  values  for 
permeability.  An  examination  of  the  data  obtained  by 
Dolch(1956)  for  samples  of  permeability  comparable  to  that 
used  in  this  study  suggested  that  the  (true)  physical 
permeability  at  infinite  mean  gas  pressure  is  about  20  to  35 
percent  lower  than  that  obtained  at  a  mean  pressure  of  one 
atmosphere.  This  variation  is  comparable  to  the  usual 
variability  in  the  measured  permeabilities;  thus,  it  was 
decided  to  not  include  such  a  correction. 
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Figure  20  Air  Flow  Velocity   Vs  Pressure  Gradient  for  S5MD 
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Figure  22  Air  Flow  Velocity  Vs  Pressure  Gradient  for  S7HD 
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4-3.2  Water  Permeability 

The  water  permeability  values  used  in  this  study  for 
silt-kaolin  mixtures  were  obtained  by  Garcia- 
Bengochea( 1978) .  He  used  a  closed  system,  falling  head 
apparatus  to  measure  the  permeability.  For  some  samples  the 
permeability  changed  with  hydraulic  gradient  (S7L0,  S7H0  AND 
S5L0) ;  thus,  Darcy's  law  was  not  obeyed  for  those  samples. 
Experimental  errors  were  given  as  the  reasons.  The 
permeability  values  assigned  to  each  of  the  samples  were 
calculated  from  an  average  of  six  to  eight  measurements  from 
the  stabilized  portion  of  the  permeability-hydraulic 
gradient  curves.  Even  though  Garcia— Bengochea  did  not 
mention  it  in  his  thesis,  particle  migration  at  the 
beginning  of  the  test  could  have  resulted  in  lower  values  of 
permeabi 1 i ty . 

For  samples  whose  compaction  variables  could  not  be 
replicated  in  this  study,  the  water  permeability  values  were 
obtained  by  calculating  the  pore  size  parameters  and  then 
using  the  prediction  equation  developed  by  him  to  calculate 
permeabi 1 i ty . 

The   water   permeability   values   of   the  natural  soils 

tested  were  obtained  by   Janssen   and   Dempsey ( 1980) .  They 

used    a    constant    head    permeability    apparatus.  The 

permeability  assigned  to  each  soil  is  an  average  from  three 
spec  imens . 
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4-4  Pore  Size  Distribution 

4-4.1  Data  Reduction 

The  pore  size  distribution  was  determined  from  the 
pressure  and  intrusion  data  using  Equation  (1-18).  The  use 
of  this  equation  requires  the  knowledge  of  the  surface 
tension  and  contact  angle  of  the  mercury.  The  surface 
tension  was  assumed  to  be  484  dynes/cm.  A  contact  angle  of 
147°  was  used  for  the  silt-kaolin  mixture  and  natural  soils 
of  high  kaolin  content.  For  the  soils  of  high 
montmor i 1 oni te  content  a  contact  angle  of  139°  was  used 
(Diamond,  1970). 

In  order  to  reduce  the  mercury  intrusion  data,  the 
computer  program  developed  by  Reed(1977)  was  used.  This 
program  includes  correction  for  compression  of  air  and 
mercury  during  intrusion.  This  program  was  modified  in 
order  to  plot  both  the  porosity  density  function  and  the 
normalized  cumulative  distribution  function  as  will  be 
discussed  in  the  next  section.  A  listing  of  the  program 
used  is  given  in  Appendix  C. 

4-4.2  Presentation  of  Data 

The  data  from  the  mercury  porosimeter  is  usually 
presented  in  the  form  of  cumulative  and  differential  pore 
size  distribution  curves.    A   detailed   discussion   on   the 
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presentation  of  data  can  be  found  in  Garcia-Bengochea( 1978) 
and  Juang(1981).  In  the  present  study  the  following  method 
was  chosen  to  present  the  data. 

a)  Cumulative  distribution  curve:  at  a  given  diameter, 
the  volume  intruded  is  normalised  to  100  percent  of 
the  total  pore  volume  and  plotted  against  the 
logarithm  of  the  diameter.  The  use  of  total  pore 
volume  to  normalize  the  distribution  is  justified  by 
the  fact  that  more  than  95  percent  of  the  pore 
volume  was  intruded  in  the  majority  of  the  samples 
tested . 

b)  Differential  distribution  curve:  even  though  the 
pore  diameter  is  a  continuous  variable,  the  mercury 
porosimeter  permits  measurements  at  discrete  points 
only.  However,  the  porosity  density  function 
obtained  from  discrete  data  is  used  in   this   study. 

The  porosity  density  function  is  obtained  from  the 
discrete  data  as  follows:  each  pressure  increment 
corresponds  to  a  decrease  in  pore  diameter  from  di_i  to  d£  . 
The  pressure  increments  are  selected  so  that  the  difference 
between  log  dj  and  log  dj.i  is  a  constant,  i.e. 

log  di  -  log  di_i  =  c  (4-1) 

di-i/di  =  10^  (4-2) 

where        d   =   pore   diameter   corresponding   to   an 

intrusion  pressure  P 
c  =  constant 
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From  Equation  (1-18),   it   is   evident   that   intrusion 
pressure  is  inversly  proportional  to  the  pore  diameter. 
Therefore , 

Pi/Pi-i  =  10^  (4-3) 

By  increasing  the  pressure  of  intrusion  according  to 
Equation  (4-3)  and  maintaining  c  constant,  logarithmically 
equal  class  width  can  be  obtained.  The  volume  of  pores 
intruded  between  diameter  dj  and  d£_i  is  calculated  from  the 
porosimeter  readings.  This  incremental  volume  is  divided  by 
the  class  width(  c  ),  normalised  to  100  percent  of  the  total 
pore  volume,  and  plotted  at  the  midpoint  of  the  class  width. 
The  midpoint  D  wi 1 1  be  obtained  by  the  following  equation, 

log  D  =  (log  di+log  di.i)/2  (4-4) 

D  =  V  d^^iTI  (4-5) 

The  porosity  density  function,  f (d) ,  can  be  obtained  by 
connecting  the  values  of  D,  f(D)  by  straight  line  segments. 
The  smaller  the  value  of  c  used  the  better  the  actual 
density  function  is  approximated.  In  this  study  a  value  of 
0.18  was  used  for  c,  Juang(1981)  found  that  the  finite 
different  approximation,  as  used  in  this  study,  compares 
well  with  the  density  function  obtained  from  direct 
differentiation  of  the  cumulative  distribution  curve. 

Garcia-Bengochea( 1978)  plotted  the  porosity  frequency 
and  suggested  that  distortion  in  the  curves  may  occur  unless 
logarithmically  equal  pressure  increments   during   intrusion 
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are  used.  In  this  study  it  was  found  that  by  plotting  the 
porosity  density  function  the  distortion,  which  could  arise 
from   unequal  pressure  increments,  can  be  kept  to  a  minimum. 

The  pore  size  parameters  were,  however,  calculated  from 
the  discrete  values  obtained  directly  from  the  porosimeter 
measurements.  This  was  done  in  light  of  the  fact  that  the 
numerical  integration  of  the  porosity  density  function  to 
get  the  poi-e  size  parameters  did  not  improve  the  accuracy 
substantially. 

4-4.3  Comparison  of  Pore  Size  Distribution  Curves 

The  pore  size  distribution  curves  are  shown  on  Figures 
24  through  26  and  in  Appendix  A  on  Figures  42  through  60. 
It  can  be  seen  that  about  two-thirds  of  the  samples  tested 
exhibit  a  bi-modal  characteristic  and  the  remaining  one- 
third  of  the  samples  show  single  modal  characteristics.  Bi- 
modal  and  single  modal  characteristic  were  found  in  both  the 
natural  soils  and  the  artificial  mixture  of  silt  and  kaolin. 

For  the  silt-kaolin  mixtures  showing  the  bi-modal 
charaterist ic ,  the  smaller  pore  mode  occured  closer  to  0.1 
or  0.2  lJ.m  and  the  larger  pore  mode  occured  between  0.5  and 
lO^m.  However,  the  majority  of  the  silt-kaolin  mixtures 
were  found  to  have  single  modal  characteristic,  and  the  pore 
mode  occured  between  0.1  and  O.Q{J.m  except  for  soil  S9MW  for 
which  it  occured  at  O.Zflm.       The  findings   for   soils   having 
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24   Pore  Size  Distribution  Curves  for  70%  Silt-30%  Kaolin 
Compacted  at  Optimum  Water  Content  for  Three  Efforts 
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Figure  25   Pore  Size  Distribution  Curves  for  50%  Silt-50%  Kaolin 
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Figure  26     Pore  Size  Distribution  Curves  for  50%  Silt-50%  Kaolin 
Compacted  at  Low  Effort 
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bi-modal  characteristic  generally  agree  with  those  of 
Garc ia-Bengochea( 1978)  whose  soils  were  all  silt— kaolin 
mixtures.  He  stated  that  all  of  his  soils  exhibited  bi- 
modal  characteristics;  however,  a  close  examination  of  his 
data  suggest  that  some  of  his  soils  did  not  exhibit  a  very 
distinct  bi— modal  characteristic. 

All  the  natural  soils  tested  in  this  research  exhibited 
bi— modal  characteristic  except  one.  The  smallest  pore  mode 
occured  between  0.4  and  1.0fJ.m,  and  the  largest  pore  mode 
occured  between  10  and  100  /zm.  For  the  soil  showing  the 
single  modal  characteristic,  the  pore  mode  occured  at  0.5/im. 
White(1980)  found  that  the  natural  soil  he  tested  exhibited 
single  modal  characteristics;  however,  he  tested  only  one 
type  of  natural  soil.  No  generalization  about  the  bi-modal 
characteristic  of  natural  soils  is  believed  appropriate 
because  even  one  of  the  natural  soils  tested  in  this  study 
exhibited  single  modal  characteristics.  This  question  is 
intriguing  and  does  deserve  further  study. 

Figure  24  shows  the  effect  of  different  compaction 
energy  levels  at  the  respective  optimum  water  contents  for 
70%  silt-30%  kaolin  mixture.  The  magnitude  of  the  peak  at 
the  larger  pore  mode  is  reduced  as  a  result  of  higher 
compaction  energy;  also,  the  magnitude  of  the  peak  at 
smaller  pore  mode  is  increased.  This  seems  to  support  the 
concept  of  the  deformable  aggregate  model.  The  smaller  pore 
mode  is  not  very  distinct   for   sample   S7L0 ;   however,   the 
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smaller  pore  mode  remains  close  to  0.1/im.  The  air 
permeability  decreases  with  increasing  energy  level,  but  the 
values  of  permeability  are  not  very  different. 

Figures  25  and  26  show  the  effect  of  water  content  at 
the  same  energy  level  for  50%  silt-50%  kaolin  mixture  at  two 
different  compaction  energy  levels.  When  the  water  content 
is  increased  the  magnitude  of  the  peak  at  pore  mode  is 
reduced  for  S5M0.  At  low  compaction  energy  level  the 
magnitude  of  the  peak  for  S5L0  is  reduced  but  shifted  to  the 
right.  This  may  be  due  to  the  sample  disturbance  as  it  was 
not  possible  to  cut  a  small  undisturbed  sample  of  S5L0 
because  of  its  higher  water  content.  The  air  permeability 
at  water  content  dry  of  optimum  is  slightly  higher  than  that 
at  optimum. 

Generally  the  results  obtained  for  silt-kaolin  mixture 
followed  the  trend  of  the  results  from  Garcia- 
Bengochea(1978) . 

4-5  Relation  between  Pore  Size  Distribution  and  Permeability 

Previous  researchers  including  Taylor(1948)  and 
Lambe(1951)  attempted  to  correlate  logarithm  of  permeability 
with  void  ratio.  Garc ia— Bengochea( 1978)  showed  that  such  a 
relationship  is  not  valid  for  compacted  fine  grained  soils. 
It  can  be  seen  from  Figure  27  that  such  a  relationship  also 
is  not  valid  for  air  permeability. 
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The  models  used  in  this  study  to  predict  the 
permeability  are  discussed  in  Section  1-1.4.  The 
relationship  between  physical  permeability  (K)  and  pore  size 
parameter(PSP)  is  generally  of  the  form 

K  =  Cg  PSP 

where  C^  is  a  shape  factor. 

The  formulas  used  to  calculate  the  pore  size  parameters 
are  shown  in  Table  1.  The  integration  was  carried  out  both 
numerically  and  from  discrete  data  using  the  method 
suggested  by  Garc ia— Bengochea( 1978) .  It  was  found  that  the 
values  of  pore  size  parameters  obtained  by  these  two  methods 
did  not  differ  significantly.  Therefore  it  is  believed  that 
as  long  as  the  logarithmic  pressure  (or  diameter)  increment 
is  kept  reasonably  small  the  pore  size  parameters  can  be 
calculated  from  discrete  data.  In  this  study  a  value  of 
0.18  was  used  for  the  exponent  c  in  the  equation  Pi/Pi_i  = 
lO'^  . 


As  can  be  seen  from  Equations  1-4  and  1-7,  the  pore 
size  parameters  are  very  sensitive  to  the  diameter  of  the 
larger  pores.  Some  of  the  larger  pores  may  be  just  cavities 
instead  of  being  continuous  pores.  Therefore,  it  was 
necessary  to  use  a  cut— off  value  for  the  larger  pore 
diameter.  A  possible  criterion  to  select  such  a  diameter  is 
to  find  a  diameter  above  which  the  contrbution  to  the  flow 
from  those  pores  is  negligible.  Since  the  soils  used  in 
this  study  are  of  different  types  it   was   not   possible   to 
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find  a  constant  cut-off  diameter  for  all  the  soils. 
Therefore,  it  was  decided  to  discard  5  to  10  percent  of  the 
pore  volume  on  the  larger  pore  sizes  in  the  calculation. 
The  pore  size  parameters  obtained  by  omitting  5  percent  of 
the  volume  did  not  differ  significantly  from  those  obtained 
by  omitting  10  percent  of  the  volume.  However,  in  this 
study  10?  of  the  pore  volume  was  omitted  from  the 
calculations  as  the  pore  size  parameters  obtained  by  this 
method  gave  a  better  correlation  with  permeability.  The 
values  of  pore  size  parameters  are  given  in  Table  7  along 
with  the  values  of  physical  permeability  of  each  of  the 
soi 1 s . 

The  pore  size  parameters  from  the  Capillary  and  the 
Marshall  model  are  plotted  against  the  physical  permeability 
obtained  from  air  permeability  tests  on  Figures  28  and  29 
respectively.  All  axes  are  on  a  logarithmic  scale.  It  can 
be  seen  that  the  relationship  is  linear.  A  linear 
regression  analysis  was  used  to  quantify  the  relationship, 
and  it  can  be  expressed  as: 

log  K  =  b  log  PSP  +  a  (4-6) 

where  a  and  b  are  linear  regression  parameters. 
Equation  (4-6)  can  be  written  as: 

K  =  10^  PSPt-  (4-7) 

According  to  the  models  used,  b  is  equal  to  1  and  10= 
is  the  shape  factor.  Table  8  presents  the  linear  regression 
equations   and  the  values  of  coefficient  of  correlation  (r^) 
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Table  7   Permeability  and  Pore  Size  Parameters 


SAMPLE 
CODE 

CAPILLARY 

MODEL 

MARSHALL 

MODEL 
xl0-2/im2 

xlO-  ^5  ni2 

S9HW 

2.28 

16.  10 

13.2 

2.25 

S9MW 

2.09 

28.70 

42.3 

3.53 

S7H0 

0.34 

0.574 

5.  17 

0.225 

S7M0 

0.27 

1.68 

7.32 

0.409 

S7L0 

0.24 

2.37 

8.82 

0.256 

S5M0 

0.021 

0.062 

0.56 

0.027 

S5MD 

0.  025 

0.059 

1.40 

0.039 

S5LD 

0.056 

0.438 

4.27 

0.  037 

S5L0 

0.  195 

1  .390 

3.20 

0.  164 

ALVIN-B 

35.9 

132.  0 

138 

2.24 

AVA-B 

1.205 

4.44 

13.6 

0.  107 

BRYCE-A 

1021 

4958 

3510 

7.80 

CISNE-B 

1661 

5250 

2560 

29.3 

HERRICK-A 

25.5 

93.0 

219 

0.017 

HOSMER-A 

12.8 

52.0 

97.8 

5.27 

HOSMER-B 

556.  1 

1460 

990 

0.654 

RICHVIEW-B 

30.5 

67.  0 

162 

0.  127 

SABLE-B 

282.2 

390 

440 

5.85 

STOY-A 

10.3 

72.0 

68.  0 

5.07 

K3=Physical  Permeability  From  Air  Permeability  Tests 
y  =Physical  Permeability  From  Water  Permeability  Tests 
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Figure  28    Physical  Permeability  (Air)    Vs    Capillary  Model 
Equivalent  Pore  Size  Parameter 
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Figure  29    Physical  Permeability  (Air)    Vs    Marshall 
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TABLE  3   LINEAR  REGRESSION  EQUATIONS  FOR  PERMEABILITY  AND 
PORE  SIZE  PARAMETERS  (PSP) 


Capillary  Model: 

K  =  (PSP)°-^2  10"^'^^ 


r^  =  0.971 


Marshall  Model: 

K  =    (PSP)O-^^    10-3-94 
r2   =   0.969 

Units: 

K  -  m^ 
PSP   -  m^ 
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for  the  models  used.  The  values  of  r^  are  close  to  one  for 
both  the  models.  This  demonstrates  the  superiority  of 
models  based  on  pore  size  distribution  to  those  based  on 
gross  volumetric  parameters.  However,  an  experimental 
matching  factor  is  necessary  to  allow  for  the  different  pore 
shapes  and  tortuosity. 

When  the  values  of  physical  permeability  obtained  from 
water  permeability  tests  were  plotted  against  the  pore  size 
parameters,  no  linear  relationship  was  obtained.  However, 
when  only  the  silt-kaolin  mixture  was  considered  a  linear 
plot  was  obtained.  This  is  probably  due  to  the  low  activity 
of  kaolin  as  will  be  discussed  in  the  next  section. 

4-6  Effect  of  Minerology  on  Permeability 

The  possible  effects  of  clay-water  interaction  on 
permeability  are  discussed  in  Sections  1-1.1  and  1-1.2.  The 
ratio  of  the  physical  permeability  obtained  by  air  and  water 
permeability  tests  reflects  the  degree  of  clay-water 
interaction.  This  ratio  for  the  soils  tested  varied  from  5 
to  1514,  and  for  the  soil  Herrick-A  it  was  12880.  For  silt- 
kaolin  mixtures  the  range  in  permeability  ratio  was  found  to 
be  between  5  and  36.  This  is  due  to  the  fact  that  the  silt- 
kaolin  mixture  has  a  low  Activity;  hence,  the  minero 1 ogica 1 
effect  on  water  during  permeation  is  very  low.  For  the  90? 
silt-10?  kaolin  mixture  the  ratio  was  between  5  and  12.  It 
should   be   noted  that  the  air  and  water  permeabilities  were 
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not  obtained  from  the  same  sample.  Thus,  the  variability  in 
the  results  and  the  particle  migration  during  water 
permeability  tests  could  have  led  to  a  difference  of  one 
order  of  magnitude  in  the  physical  permeabilities  obtained 
from  air  and  water  permeability  tests.  On  the  other  hand, 
for  the  natural  soils  tested  the  permeability  ratio  varied 
from  13  to  12880.  These  soils  have  Activities  ranged  from 
0.13  to  2.2.  Also,  X-ray  diffraction  analysis  showed  that 
some  of  these  soils  have  traces  of  montmor i 1 1 oni te  which  is 
a  very  active  mineral. 

It  has  been  believed  that  the  effect  of  adsorbed  layer 
is  to  reduce  the  sizes  of  the  pores  available  for 
permeation,  thus  reducing  the  permeability.  To  check  this, 
pore  size  parameters  were  calculated  using  pore  sizes 
reduced  by  an  assumed  thickness  for  adsorbed  layer,  and 
these  parameters  were  substituted  into  the  predictive 
equation  already  obtained  to  calculate  permeability.  The 
permeability  thus  obtained  was  much  higher  than  that 
measured.  As  discussed  in  Section  1-1.2,  earlier 
researchers  also  found  that  the  adsorbed  layer  does  not 
affect  the  permeability  significantly.  Therefore,  it  was 
concluded  that  the  adsorbed  layer  would  not  have  caused  such 
a  large  difference  in  permeability  with  water  and  air.  The 
e Ictro-kinetic  coupling  was  also  found  to  be  insignificant 
by  various  researchers  including  01sen(1962). 
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As  discussed  in  Section  1-1.2,  Mesri  and  01son(1971) 
also  found  large  differences  in  permeability  of  clays  to 
polar  and  non-polar  fluids.  The  mechanism  suggested  by  them 
can  be  extended  to  the  soils  tested  in  this  study  as 
follows:  The  soil  particles  form  stable  aggregates  and  the 
primary  fluid  flow  takes  place  in  the  larger  voids  between 
the  aggregates.  In  the  presence  of  water,  which  is  a  polar 
fluid,  the  aggregates  break  down  and  a  dispersive  structure 
of  low  permeability  results.  The  degree  of  dispersion 
depends  on  the  physico-chemical  properties  of  soil  minerals. 
Since  the  increase  in  activity  reflects  the  presence  of 
highly  active  minerals,  a  correlation  between  permeability 
ratio  and  activity  was  attempted.  The  values  of  activity 
and  permeability  ratio  of  the  soi.s  tested  are  given  in 
Table  9  and  Figure  30  shows  the  plot  of  activity  vs 
logarithm  of  permeability  ratio.  Even  though  considerable 
scatter  is  found  in  the  plot,  it  is  useful  in  predicting  the 
water  permeability  of  active  soils  from  pore  size 
distribution  measurements,  because  all  the  models  proposed 
so  far  assume  that  the  pore  size  distribution  does  not 
change  during  flow. 

4-7  Soil  Moisture  Characteristic  Curves 


In  this  study,  the  mercury  intrusion  technique  was  used 
to  obtain  the  soi 1 -moisture  characteristic  curve.  The 
theory   behind   this   method   is  described  in  Section  1-2.3. 


Table  9   Activity  and  Permeability  Rati( 
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SAMPLE 
CODE 

ACTIVITY 

PERMEABILITY 
RATIO 
Ka/Kw 

S9HW 

NA 

5.SS 

S9MW 

NA 

11.81 

S7H0 

0.33 

22.9 

S7M0 

0.33 

17.9 

S7L0 

0.33 

34.4 

S5M0 

0.29 

20,7 

S5MD 

0.29 

36.0 

S5LD 

0.29 

11.5 

S5L0 

0.29 

19.5 

ALVIN-B 

0.  13 

61.6 

AVA-B 

0.48 

127.0 

BRYCE-A 

2.2 

450 

CISNE-B 

0.69 

87.3 

HERRICK-A 

0.92 

12880 

HOSMER-B 

0.73 

1514 

RICHVIEW-B 

0.38 

1275 

SABLE-B 

0.87 

75.2 

STOY-A 

1.3 

13.4 
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Figure  30    Activity  Vs  Permeability  Ratio 
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The  predicted  curves  were  compared  with  the  ones  obtained 
from  experimental  method  by  Janssen  and  Dempsey ( 1980)  for 
the  same  soils.  They  used  the  Tempe  pressure  cell  and 
pressure  plate  apparatus  to  obtain  the  soi 1— moisture 
characteristic  curve.  The  highest  value  of  suction  measured 
was  15  bar (1500  KN/M^).  Figures  31  through  41  show  the 
soi 1— moisture  characteristic  curves  obtained  by  experimental 
method  and  those  predicted  from  pore  size  distribution 
measurements . 

It  can  be  seen  from  Figures  31  through  36  that  the 
soil-moisture  characteristic  curve  is  very  well  predicted  up 
to  a  suction  oC  3000  cm  of  water  whereas  for  soils  ALVIN-B, 
HERRICK-A  and  STOY-A  the  prediction  is  good  only  up  to  1000 
cm  of  water.  However,  considerable  difference  between  the 
experimental  and  predicted  curves  were  found  for  soils  shown 
in  Figures  37  through  39.  An  examination  of  the  values  of 
the  compacted  variables  revealed  that  considerable 
difference  in  dry  density  and  water  contents  existed  between 
the  samples  used  in  the  mercury  intrusion  method  and  the 
conventional  pressure  plate  test.  Also  the  experimental 
curves  were  not  replicated.  Thus  the  variability  in  results 
could  have  caused  these  differences. 

The  Figures  40  through  41  show  that  above  a  suction  of 
3000  cm  of  water  the  experimental  and  the  predicted  curves 
diverge  considerably.  This  can  be  expected  because  water 
retention   in  the  higher  suction  range  is  controlled  more  by 
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Figure  32   Soil  Moisture  Characteristic  Curves  for  Cisne-B 
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Figure  33   Soil  Moisture  Characteristic  Curves  for  Hosmer-A 
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Figure  34  Soil  Moisture  Characteristic  Curves  for  Richvlew-B 
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Figure  35   Soil  Moisture  Characteristic  Curves  for  Sable-B 
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Figure  36    Soil  Moisture  Characteristic  Curves  for  Alvin-B 
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Figure  37    Soil  Moisture  Characteristic  Curves  for  Herrick-A 
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Figure  38   Soil  Moisture  Characteristic  Curves  for  Stoy-A 


119 


lOOOOQ 


lOOOQ 


CC 
2: 


O 


O 


O 

ID 
CO 


lOOQ 


100 


10    .. 


SIMPLE  CODE        BRYa-fl 

HflTER  aSTTENTtCS'T.)  =     2H.t 
DRY  DENSITY  =         87.2   LB/OI.FT 

♦    EXPERlMENTfL  METHOD 

(JflNSSEN  AND  DEHPSEY.  1980) 


FROM  PSD  MEASUREMENTS 


1      I     ,     ■     ■     . 


I     I     t     I     I     I     I     I     t — I     I     I     I     I     I     I     I — I     I     I 


CI 
0       CD 

® 


|] 


-w 


m 

n 

Q 

o    m 


a 


10  15  20  25 

WATER  CONTENT  iZ 


30 


35 


Figure  39    Soil  Moisture  Characteristic  Curves  for  Bryce-A 
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Figure  40    Soil  Moisture  Characteristic  Curves  for  Catlin-B 
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adsorption  and  less  by  structure.  The  mercury  intrusion 
method  does  not  consider  such  physico-chemical  effects. 
However,  the  ability  to  predict  suction  values  as  large  as 
3000  cm  of  water  is  adequate  for  practical  purposes. 
Moreover,  the  time  required  for  the  mercury  intrusion  method 
is  minutes  as  compared  to  many  days  for  conventional  methods 
to  determine  soi 1 -moisture  characteristic  curves. 
Therefore,  this  study  demonstrates  that  the  equilibrium 
water  content  can  be  predicted  and  the  basis  for  prediction 
is  readily  obtained. 

4-8  Frost  Heave 

At  the  beginning  of  this  study,  it  was  hoped  to 
correlate  the  frost  heaving  rate  of  compacted  silt-kaolin 
mixtures  with  effective  pore  size  distribution.  The  frost 
heaving  rates  for  such  soils  were  obtained  by  Reed(1977). 
However,  as  disscussed  in  Section  4-6,  it  was  found  that  the 
adsorbed  layer  does  not  affect  the  flow  of  water 
considerably  for  silty  clays.  Therefore,  the  predictive 
equation  for  silty  clays  as  determined  by  Reed(1977) 
requires  no  modifications. 

For  clays  in  which  permeability  is  affected  by  the 
clay— water  interaction  further  research  is  required  to 
establish  relationship  among  pore  size  distribution,  frost 
heaving  rate,  and  physico-chemical  variables. 
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5  -  CONCLUSIONS  AND  RECOMMENDATIONS 

5—1  Cone lusions 

This  study  concerned  artificial  mixtures  of  silt  and 
kaolin  as  well  as  samples  of  naturally  existing  soils. 
Recognizing  that  small  specimens  were  tested  in  the 
laboratory,  the  following  conclusions  are  offered: 

1)  Bi-modal  and  single  modal  pore  size  distributions 
occur,  but  no  association  could  be  made  which  would 
either  explain  or  predict  the  modality. 

2)  The  modality  of  the  pore-size  distribution  curves 
appears  to  be  real.  Using  different  numerical 
procedures  for  curve  plotting  and  parameter 
determination  did  not  markedly  improve  correla- 
tions. A  differential  pore-size  distribution  curve 
in  the  form  of  the  porosity  density  function  is 
believed  clearer  for  the  observer. 

3)  There  appears  to  be  no  single  actual  "effective 
diameter"  which  characterizes  the  pore-size 
distribution  sufficiently  for  correlation  purposes. 
A  combination  of  pore-size  distribution 
characteristics  appears  necessary  to  adequately 
describe  the  distribution. 

4)  Changes  in  the  compaction  variables   with   which   a 
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given  soil  is  prepared  produce  discernible  and 
consistent  changes  in  the  pore  size  distribution 
characteristics.  This  confirms  the  results  of  the 
previous  researchers  Garcia-Bengochea( 1978)  and 
Juang(1981). 

5)  Darcy's  law  is  obeyed  for  flow  of  gases  through 
porous  media.  The  flow  velocity  varies  linearly 
with  the  pressure  gradient.  However,  a  small 
threshold  gradient  is  found  to  exist,  but  it  is 
believed   to   be  the  result  of  experimental  errors. 

6)  Air  permeability  can  be  predicted  from  pore-size 
distribution  parameters.  The  models  for  the 
prediction  do  not,  however,  also  predict  water 
permeability  with  any  accuracy.  A  correlation 
between  air  and  water  permeability  was 
accorap 1 ished. 

7)  Predictive  equations  for  permeability  based  upon 
grain-size  parameters  are  unsuitable  for  fine 
grained  soils. 

8)  The  inability  of  models  for  air  permeability  to 
predict  water  permeability  is  a  manifestation  of 
the  effect  of  soil  minerology,  especially  the  clay 
content  and  composition.  The  larger  the  Activity 
of  the  soil,  the  larger  appears  to  be  the 
difference  between  air  and  water  permeability. 

9)  It  is  believed  that  the  effect  of  minerology  upon 
permeability    is   through   the   development   of   a 
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dispersed  structure  during  permeation. 

10)  The  soil  moisture  characteristic  curve  (ie. 
equilibrium  soil  water  content  as  a  function  of 
soil  suction)  can  be  readily  predicted  from  the 
pore  size  distribution  data.  The  time  required  to 
obtain  the  data  is  minutes  as  compared  to  many  days 
for  suction  data.  The  prediction  thus  made  appears 
to  be  adequate  for  suctions  as  large  as  about  25 
meters  of  water. 

11)  The  in-situ  equilibrium  water  content  of  a  compacted 
soil  can  be  predicted  provided  that  an  estimate  is 
possible  of  the  location  of  the  free  water  level  in 
the  field  installation. 

13)  Pore-size  data  appears  to  be  a  sensitive  vehicle 
for  monitoring  behavior  effects  that  are  controlled 
by  soil  fabric. 
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5-2  Recommendations  for  Future  Research 

Further  research  is  needed  in  the  following  areas: 

1)  Obtain  a  correlation  among  permeability  ratio,  pore 
size  parameters,  and  physico-chemical  variables 
such  as  ion  exchange  capacity. 

2)  In  this  study  a  linear  relationship  was  obtained 
between  air  flow  rate  and  pressure  gradient. 
Professor  G. Bl ight ( 1981-pr ivate  communication) 
suggested  that  the  fixed  head  apparatus  tends  to 
give  a  non— linear  relationship.  This  difference 
should  be  examined. 

3)  Influence  of  pore  size  distribution  on  strength  and 
compressibility  requires  attention. 

4)  Effect  of  ion-exchange  capacity  on  the  soil- 
moisture  characteristic  curve  at  higher  suction 
values  should  be  examined. 

5)  The  minerology  effect  on  frost  heave  rate  should  be 
studied  for  highly  active  clays. 
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APPENDIX  A 


Pore  Size  Distribution  Curves 
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Figure  42    Pore  Size  Distribution  Curves  for  Sample  S9HW 
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Figure  44    Pore  Size  Distribution  Curves  for  Sample  S7H0 
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Figure  45    Pore  Size  Distribution  Curves  for  Sample  S7M0 
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Figure  46    Pore  Size  Distribution  Curves  for  Sample  S7L0 
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Figure  48    Pore  Size  Distribution  Curves  for  Sample  SSMO 
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Figure  49    Pore  Size  Distribution  Curves  for  Sample  S5LD 
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Figure  52    Pore  Size  Distribution  Curves  for  Sample   Ava-B 
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Figure  53    Pore  Size  Distribution  Curves  for  Sample  Bryce-A 
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Figure  54    Pore  Size  Distribution  Curves  for  Sample  Cisne-B 
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Figure  55    Pore  Size  Distribution  Curves  for  Sample  Herrick-A 
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Figure  56    Pore  Size  Distribution  Curves  for  Sample  Hosmer-A 
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Figure  57    Pore  Size  Distribution  Curves  for  Sample  Hosmer-B 
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Figure  58    Pore  Size  Distribution  Curves  for  Sample  Richview-B 
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Figure  59    Pore  Size  Distribution  Curves  for  Sample  Sable-B 
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Figure  60   Pore  Size  Distribution  Curves  for  Sample  Stoy-A 
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APPENDIX  -  B 

Derivation  of  Permeameter  Equation  (from  Dolch,  1956) 
Figure  6  1  shows  the  schematic  diagram  of  the  air  permeability 

apparatus  used  in  this  study. 

The  mass  of  air,  m,  on  the  low  pressure  side  of  the  sample  at  any 

time  during  the  decay  process  is 

m  =  P^CVh  -  h  ^)pi  (Bl) 

2 

where    Ph  -  Pressure  on  the  low  pressure  side  of  the  sample 

V,  -  Volume  of  the  apparatus  on  the  low  pressure  side 

of  the  sample  when  the  differential  manometer  is 

at  equilibrium 

h  -  Height  difference  on  the  differential  manometer 

A  -  Internal  cross  sectional  area  of  the  differential 
t 

manometer  tube 
p   -  Density  of  the  air  at  unit  pressure,  i.e.  p  =  m/RT 
where     m  =  Molecular  weight 
R  -  Gas  constant 
T  -  Absolute  temperature 
The  time  rate  of  change  of  m  is  the  mass  flow  rate  assuming  steady 
state  conditions.   These  are  not,  of  course,  steady  state  conditions. 
Kamack  (1954),  however,  stated  for  a  similar  set  of  circumstances 
"...a  negligible  error  is  involved  in  assuming  that  the  air  instantly 
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Figure  6  I     Schematic  Diagram  of  the  Air  Permeability  Apparatus 
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adjusts  itself  to  the  equilibrium  flow  rate  corresponding  to  the 
pressure  drop.   This  situation  is  usually  referred  to  as  a  quasi- 
stationary  state".   If,  then,  steady  state  conditions  are  assumed  and 
equation  (Bl)  is  differentiated  with  respect  to  t,  the  result  is 


dm 
dt 


=  PI 


dPk  _  :!l  ( !!^  +  ^hdh  )] 

_dt    2     dt      dt  J 


The  approach  velocity  at  mean  pressure,  u^,  is  given  by 

Vl^m^dt^ 
where     A  =  Approach  area  of  the  sample 

P  =  Mean  pressure  across  the  sample 


(B2) 


(B3) 


From  B2  and  B3 


^a  =  ap: 


'   dPh   At   hdPh  Phdh  n 
_^h  dT"  "  2"  ■  "dF  ^   IF  ^J 


(B4) 


The    physical   permeability,   K,    for   the   case   of   gaseous    flow   is 


defined   by  Muskat ,   via  Darcy's    law,    as 

K  AP 


y  L 


(B5) 


where 


AP  =  head  loss  across  the  sample  in  pressure  units 
y  =  viscosity  of  the  gas 
L  =  macro  length  of  the  sample 
substituting  B4  in  B5  yields, 


^   AP  As  Pm 
In  SI  units, 

y  =  KN.s/m^ 
AP,Ph,Pni  =  KN/m2 

V,  =  m^ 
h 


dPv 


h"t,  "^h,   Ph  ""t  dh 
(Vh  -  — ^^dF"^  -    2  dT 


(B6) 
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L,h   =  m^ 


^s'^t    =   "^ 


t    =   sec 


K  =  m^ 
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APPENDIX  C 


Pore  Size  Distribution  Computer  Programs 
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PROGRAM  PORE  CIMPUT, OUTPUT. TAPE5=iriPUT,TAPEG=0UTPUT, PLOT) 
C 

C     THIS  PROGRAM  CALCULATES  AND  TABULATES  PORE  SIZE  DISTRIBUTION 

C     PARAMETERS,  ALONG  WITH  PLOTTING  PORE  DIAMETER  UERSUS  CUMULATIUE 

C     POROSITY  AND  POROSITY  DENSITY  FUNCTION.  THIS  PROGRAM  WAS 

C     ORIGINALLY  WRITTEN  BY  REED(1977)  AND  MODIFIED  BY  PRAPAHARAN(1977) 

C 

INTEGER  TITLEO) 

REAL  CMSP(50),CUPORO(50) 

REAL  S(50),P(50).PS(50),D(50),SP(50).CUSP(50) 
C 

A=0.0 

READ.NSET 

CALL  PLOTS 

DO  22        L=1,NSET 

NL=0 

B=0.0 

READ.NDATA 
C 

DO  21   J=1,NDATA 

READ(5. 1)  TITLE 

1  FORNATOAIO) 
WRITECG.S)  TITLE 

2  F0RMAT(1H1///5X,?!SAMPLE  CODE:   ?i,3A10///) 
C 

C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  BE  ABLE  TO  CALCULATE  RESULTS  FOR 

C     MORE  THAN  ONE  SET  OF  DATA  IN  A  SINGLE  RUN.   NDATA  IS  THE  NUMBER 

C     OF  SETS  OF  DATA  TO  BE  CALCULATED.   NDATA  IS  AN  INTEGER  UALUE.  AND 

C     IS  THE  FIRST  DATA  CARD  FOR  THE  SET  OF  DATA. 

C     DETERMINATION  OF  PORE  SIZE  DISTRIBUTION 

C 

C 

C 

READ, US, DS,DM, UP, WSP, W5PM, PE, PF, SR 

PRINT, US, DS, DM, UP, WSP, WSPM, PE, PF, SR 
C 

C     WS=WEIGHT  OF  SAMPLE,  DS=SPECIFIC  GRAUITY  OF  SOLIDS,  DM=DENSITV  OF 
C     MERCURY,  UP=UOLUME  OF  PENETROMETER,  WSP=UEIGHT  OF  SAMPLE  AND  PEN- 
C     ETROMETER,  WSPM=WEIGHT  OF  SAMPLE  ,  PENETROMETER  AND  MERCURY, 
C     PE  =EUACUATING  PRESSURE,  PF=F1LLING  PRESSURE, IN  MM  OF  MERCURY, 
C     SR  =STEM  READING  AT  FILLING  PRESSURE.   THE  ABOUE  UALUES  ARE  ALL  ON 
C     THE  SAME  DATA  CARD 
C 

U5=WS/DS 

UM=WSPM  -  WSP 

UM=WM/DM 

USA=UP-UM*( 1 . +PE/PF)-SR 

UUO=USA-US 

UOIDR=UUO/US 

POROS=UUO/USA 

PU=PE*UM 

UAI=PU/PF 
C 

WRITE(B,3)  US,UUO,UOIDR,POROS 

3  FORMAT(///T10,?^UDLUME  OF  SOLIDS  =^,F8.G/'/ 

1  Tll,i^UDLUME  OF  UOIDS  =?i,FS.E// 

2  TlG,?iUOID  RATIO  =?i,F8.G// 

3  T18,;iP0R0SITY  =?i,F8.G//) 
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READ.N,M»SI 
C 

C     M=INTEGER  UALUE  FOR  THE  NUMBER  OF  LOW  PRESSURE  REftDINCS  MADE,  NOT 
C     COUNTING  THE  FILLING  PRESSURE  READING.   M=INTEGER  UALUE  FOR  THE 
C     NUMBER  OF  HIGH  PRESSURE  READINGS  MADE,  NOT  COUNTING  THE  INITIAL 
C     POROSIMETER  READING.   SI=INTRUSION  READING  AT  INITIATION  PRESSURE 
C     IN  ML.   THE  THREE  UALUES  ARE  ALL  ON  ONE  DATA  CARD. 
C 

C     LOU  PRESSURE  CALCULATIONS 
C 

URITE(G.4) 

4  F0RMAT(1H1//,3?X,?!L0U  PRESSURE  CALCULATIONS?!///) 
WRITECG.S) 

5  FORMAT (//T5,?^MM  OF  HG^,  TEO.^PRESSURE^.TSSrPiUOLUME  0F^,J5Q,^S^E\^^f 

1  TG5.  J^INTRUSION;^,  T80,  ^^INCREMENTAL?!.  T95.  ?:CUMULATIUE^,  TllO, 

2  =!DIAMETER?:/T38.?!AIR?i, 

3  T43,?:READING?:.T81,^POROSITY;ii,TSG,?:POROSITY;^,T110,p!(MICRON)?!//) 


CUSP(l)  =  0. 


DO  7   1=2, N+1 
READ,PS(I),S(I) 
C 

C     P(I)  AND  S(I)  ARE  THE  PRESSURE  AND  INTRUSION  UALUES  FOR  EACH 
C     READING,  P(I)  IS  GIUEN  IN  MM  OF  MERCURY  FOR  THE  LOU  PRESSURE 
C     INTRUSION  AND  PSI  FOR  THE  HIGH  PRESSURE.   SCI)  IS  IN  ML.   THE  TWO 
C     UALUES  ARE  TYPED  ON  ONE  DATA  CARD,  WITH  A  CARD  FOR  EACH  READING. 
C 

P(I)=PSCI)»0. 01334 

UA=PU/PS(I) 

CUA=UAI-UA 

CSR=S(I)-SR 

SN=CSR-CUA 

IF  (SN.LT.O.)  SN  =  0.0 

SP(I)  =  SN/USA 

CUSP(I)  =  CUSP(I-l)  +SP(I) 

CUPORO(I)  =  CUSP(I)/POROS 

D(I)=((-4.)*484.»C0S(2.5G5G)».145)/P(I) 

IF  (CUPORO(I).LE.O.IO)  DM1  =  D(I) 
C 

C     P(I)=INTRUSION  PRESSURE  IN  MM  OF  MERCURY,  PS=INTRUSION  PRESSURE  IN 
C     PSI,  UA=UOLUME  OF  AIR  IN  ML,  SCI)=STEM  READING,  SN=CORRECTED 
C     INTRUSION  FOR  EACH  INCREMENT  IN  ML,  SG=CORRETED  INTRUSION  /GRAM 
C     FOR  EACH  INCREMENT,  CU(I)=CUMULATIUE  CORRECTED  INTRUSION/GRAM, 
C     SP(I)=INCRENENTAL  POROSITY,  CUSPC I)=CUMULATIUE  POROSITY, 
C     D(I)=LINITING  PORE  SIZE  DIAMETER 
C 

WRITECG.G)  PS(I),P(I),UA,S(I),SN,SP(I),CUSP(I),D(I) 
G     FORMAT(4X,F7.1,T21,F5.2,T3G,FG.5,T50,F5.4,TG7,FG.5,T83,FG.5, 
1       T97,FG.5,T111,FG.E) 

UAI=UA 

7  SR=S(I) 
C 

C     HIGH  PRESSURE  CALCULATIONS 
C 

WRITE(G,8) 

8  F0RMAT(//,3GX,;^HIGH  PRESSURE  CALCULATIONS?^//) 
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WRITE(G,3) 

FORriAT(/,T5.?iPRESSURE;^.T20.?*STEM?i,T35.?:MER.    COMP.?!.T50. /INTRUSION?!, 

1  TG5.?:iNCREMEHTAL?=,T80,?iCUMULftTIUE/,T95,':DIPiriETER?:/T7,?:(PSI)/,T19, 

2  ?:READING?:,  T35.  ^iCORRECTION?!,  TS7,  ?:POROSITY?i,  T82,  ?:POROSITYi^,  T95. 

3  :»!( MICRON)?^/) 

HCI  =  11.«.00H*UM/15000. 


DO  11   I=N+2.M+N+1 

READ, PCI), S(I) 

P(I)=P(I)+11. 
C 

C     HC  IS  CORRECTION  FOR  COMPRESSION  OF  MERCURY 
C 

HC=P(I)».0011*UM/15000. 

CHC=HC-HCI 

CSR=S(I)-SI 

SN=CSR-CHC 

SP(I)  =  SH/USA 

CUSP(I)  =  CUSP(I-l)  +  SP(I) 

CUPORO(I)  =  CUSP(I)/POROS 
C 

C     D(I)  PORE  DIAMETER  CALCULATIONS  BASED  ON  A  SURFACE  TENSION  UALUE 
C     OF  484  DYNES  /  CM  AND  A  CONTACT  ANGLE  OF  147  DEGREES. 
C 

DCI)=((-4.)*484.*C0S(2.5G5S)».145)/P(I) 

IF  (CUPORO(I).LE.O.IO)  DM1  =  D(I) 
C 

C     P(n=INTRUSION  PRESSURE  IN  PSI,  S(I)=STE!1  READING,  HC=MERCURY 
C     CORRECTION.  CHC=CHANGE  IN  MERCURY  CORRECTION  IN  ML,  SN=CORRECTED 
C     INTRUSION  IN  ML,  SGH=CORRECTED  INTRUSION  /GRAM,  CU( I )=CUMULATIUE 
C     CORRECTED  INTRUSION/GRAM,   SP(I )=INCREMENTAL  POROSITY,  CUSP(I)= 
C     CUMULATIUE  POROSITY) 
C 

WRITECG, 10)  P(I).S(I), CHC, SN, SP( I ) , CUSP( I ) , DC  I ) 

10  FORMAT(4X,F7,1,T20,F5.4,T3G,FG.5,T51,FG.5,TG8,FG.5,TS3.FB.5, 
1       T9S,FG.3) 

C 

HCI=HC 

11  91=9(1) 
C 

D(1)=-4.*484.*C05(2.5G5G)*0. 145/(0. 01934»PF) 

C 

DO  12   I  =  2,N+M+1 
CM9P(I)  =  9P(I) 

12  CONTINUE 
C 

C     BELOU  -  CALCULATION  OF  MEAN  PORE  DIAMETER  AND  SECOND  MOMENT  ABOUT 

C     THE  ORIGIN.   10  PERCENT  OF  THE  PORE  UOLUME  ON  THE  LARGER  PORE 

C     SIDE  IS  IGNORED 

C 

C 

SUMl  =0.0 

SUM2  =0.0 

DO  13   I=2,N+M+1 

DM  =  9QRT(D(I)»D(I-1)) 

IF(DM.GT.DMl)   DM=0.0 

gUMl  =  SUNl  +  SP(I)»DM 
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13  SUM2  =  SUM2   ■••  SP(I)»DM**E 
DMEA1=   SUni/POROS 
DriEAE=   SUME/PORQS 

WRITE    (G.14) 

14  FORMAT(1H1,5(/),10X,?:THE  STATISTICAL  MEASURES  BELOW  MEGLECT  PORES?: 
A,/,10X,?iMHEN  THE  CUMULATIUE  POROSITY  IS  LESS  THAIS  10  PERCENT?:///) 

URITE(G,15)   DMEAl,DNEA2,SUri2 

15  FORMAT(/T10,?:ADJUSTED  MEAN  PORE  SI2E=?:,E10.4// 
1T10,?:ADJUSTED  SECOND  MOMENT  ABOUT  THE  ORIGIN=?:,E10.4// 
2T10.?:CAPILLARY  PARAMETER  (IN  MICR0N»»2)  =  ?:,E10.4) 

C 

C     BELOW  -  CALCULATION  OF  MARSHALL  MODEL  PORE  DIAMETER 
SIGMA=0.0 
DO  17   K=2,N+M+1 
DM=D(K) 
IF  (D(K).GT.DMl)  CMSP(K)  =  0.0 
DO  IG   I=2,N+M+1 

IF(D(I1.LT.D(K))   DM=D(I) 
IF(DM.GT.DNl)     DM=0.0 
IF(D(I).GT.DM1)  CMSP(I)=0.0 
SIGMA=SIGMA+CMSP(n«CMSP(K)*DM»«2 
IB     CONTINUE 

17  CONTINUE 
WRITECG.IS)  SIGMA 

18  FORMAT(//T10,?:MARSHALL  PORE  PARAMETER  (IN  MICR0N»»2)  =?:»E10.4/) 
C 

WRITE(B,19)   DM1 

19  F0RMAT(///9X.?:CUT  OFF  DIAMETER  =?:,F7.2//) 
C 

DO  20   I=2,N+M+1 

20  D(I)=ALOG10(D(I)) 

C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  CONUERT  THE  PORE  SIZE  DIAMETER 

C     TO  THE  LOG  OF  THE  DIAMETER 

C 

CALL  PLOT(0,0.3,-3) 

CALL  JPL0T2 ( CUSP. D, M. N. A» NL. B, POROS. J. TITLE ) 

CALL  PLOT(0,-0.3,-3) 
C 

CALL  PLOTCO, 1.2.-3) 

CALL  JPLOTl ( D, SP. M. N, A, B. PF» NL, J» TITLEr POROS) 

CALL  PLOTCO. -1.2,-3) 
C 

NL=NL+1 

B=B+0.33 

21  CONTINUE 
A=A+10.0 

22  CONTINUE 

CALL  PLOT(0. 0,0. 0.999) 
STOP 
END 
C 

c 

SUBROUTINE  JPLOTl (D, SP, M, N, A, B. PF, ML. J, TITLE. POROS) 
DIMENSION  SP(50).D(50),D1(50) 
INTEGER  TITLE(3) 
CALL  FACTOR  (1.0) 
IF(J.GT.l)  GO  TO  7 
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CALL  SYMBOL  (A+1.85, .5, .7/G.,33HLIMITIMG  PORE  DIAMETER  IM  MICRONS. 
10.0,33) 

CALL  SYMBOL(A+1.,0.75,0.7/S.. 
F50H.01      0.1       1.0      10.0     100.0    1000.0.0.50) 
C     THE  ABOUE  PLOTS  AND  LABELS  THE  X  AXIS 
DO  1   1=2,5 
CI  =  I 
CALL  SYMBOL  (A+CI. 1 .0, 0. 1, 13, O.O.-l ) 

1  COMTINUE 

C     THE  ABOUE  PLACES  TICK  MARKS  ALONG  THE  X  AXIS 

CALL  SYMBOL  ( A+0.5, 1 . 15, .7/G. ,29HP0R0SITY  DENSITY  FUNCTION  (/). 
1  90.0.29) 

CALL  SYMBOL  (A+0.B05. 1. 0, .7/G. , IHO, 0.0. 1 ) 

CALL  SYMBOL  (A+0.G05,2. 0, .7/G. ,2H50, 0.0,2) 

CALL  SYMBOL  (A+0.GO5,3.0, .7/G. , 3H100, 0.0, 3) 

CALL  SYMBOL  (A+0.G05,3.9, .7/G. , 3H150. 0.0,3) 
C     THE  ABOUE  PLOTS  AND  LABELS  THE  Y  AXIS 

DO  2   1=1.4 

CI  =  I 

CALL  SYMBOL  (A+1. 0. CI. 0. 1. 15, O.O.-l ) 

2  CONTINUE 

C     THE  AEDUE  PLACES  TICK  MARKS  ALONG  THE  Y  AXIS 

CALL  SYMBOL  ( A+4.40, 3.80. .OSS.SHLEGEND,  0. 0.  G) 
C     THE  ABOUE  WRITES  LEGEND  ON  THE  PLOT 

CALL  PLOT  (A+1. 0,4. 0.3) 


CALL 

PLOT 

(A+G.0.4.0.2) 

CALL 

PLOT 

(A+G.0.2.5.2) 

CALL 

PLOT 

(A+-4.0.3.0,3) 

CALL 

PLOT 

(A+1. 0,3. 0.2) 

CALL 

PLOT 

(A+1. 0.2. 0.3) 

CALL 

PLOT 

(A+G.0.2.0,2) 

CALL 

PLOT 

(A+G. 0,1.0,3) 

CALL 

PLOT 

(A+1. 0,1.0,2) 

CALL 

PLOT 

(A+1. 0,4. 0,2) 

CALL 

PLOT 

(A+2.0,4.0,3) 

CALL 

PLOT 

(A+2. 0,1.0,2) 

CALL 

PLOT 

(A+3. 0,1.0,3) 

CALL 

PLOT 

(A+3.0,4.0,2) 

CALL 

PLOT 

(A+4.0,4.0,3) 

CALL 

PLOT 

(A+4. 0,1.0.2) 

CALL 

PLOT 

(A+5. 0.1.0,2) 

CALL 

PLOT 

(A+5.0.2.5.2) 

CALL 

PLOT 

(A+G. 0.2. 5, 3) 

CALL 

PLOT 

(A+G. 0,1.0,2) 

CALL 

PLOT 

(A+1. 0,1.0,2) 

30  4 

IOU= 

=  1.5 

CI0U= 

=  IOU 

DO  3 

10=2.9 

CIOA= 

=  10 

CIO=f 

^LOGIO(CIOA) 

CALL 

SYMBOL  (A+CIOU+CIO.1.0. 

.05.13.0 

3 

CONTINUE 

4 

CONTINUE 

CALL 

PLOT 

(A+1. 0.1.0,3) 

DO  G 

10= 

1,3 

CIO= 

10 

DO  5 

IOU= 

=  1,4 

CIOU: 

=  IOU 

-1) 
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CI0U=CI0U/-5. 

CALL  SYMBOL  (A+1.0,CIO+CIOU, .05, 15, 0.0,-1) 
5     CONTINUE 
G     CONTINUE 

7  CONTINUE 

CALL  PLOT  (A+1.0,4.0,3) 

D(1)=ALOG10 (-4. »484.«COS(2.5E5B)«0. 145/(0. 01334»PF)) 

DO  8   I=2.N+M+1 

D1(I-1)=D(I-1)-D(I) 

D(I-1)=(D(I)+D(I-1))/E 

8  SP(I-1)=SP(I)/P0R0S/D1(I-1) 
D(n+N+l)=-3.0-A 
D(M+N+2)=1.0 
SP(N+N+l)=-0.5 
SP(M+N+2)=0.5 

CALL  SYMBOL  (ft+4. 10. 5.53-B, .07,NL, 0.0,-1) 

CALL  SYMBOL  ( A+4.30, 3.55-B, . 080, UHSAMPLE  CODE,0.0.11) 

CALL  SYMBOL  ( A+5. 20, 3.55-B, . 080, TITLE, 0.0,30) 
C 

C     THE  ABOUE  LISTS  THE  WATER  CONTENT  AND  PERMEABILITY  FOR  THE  PORE 
C     SIZE  CURUE 
C 

CALL  LINE(D,SP,M+N,1,1.NL) 
C 

C     THE  ABOUE  PLOTS  THE  ACTUAL  PORE  SIZE  CURUE 
C 

C     B  IS  FOR  THE  PLACEMENT  OF  THE  PERMEABILITY  AND  WATER  CONTENT 
C     UALUES. 

RETURN 

END 
C 

C      ••»««•««»««««»««»««««•«««»«•#«•»»«  j,>f  «««•«»•«••»•«••»«•»•»•*«»•« 
C 

SUBROUTINE  JPL0T2(CUSP, D, M, N, A, NL, B, POROS, J, TITLE) 

INTEGER  TITLE(3) 

REAL  D(50),CUSP(50) 

CALL  FACTOR  (1.0) 

IF(J.GT.l)   GO  TO  3 

CALL  SYMB0L(A+1.85,5.5,.7/B.,33HLIMITING  PORE  DIAMETER  IN  MICRONS, 
GO. 0,33) 

CALL  SYMBOL(A+1.,5.75,0.7/B., 
F50H.01      0.1       1.0      10.0     100.0    1000,0.0,50) 
C     THE  ABOUE  PLOTS  AND  LABELS  THE  X  AXIS 

CALL  SYMBOL  (A+2. 0, G. 0, 0. i, 3, 0.0,-1 ) 

CALL  SYMBOL  (A+3. 0, G. 0, 0. 1, 3, 0.0,-1 ) 

CALL  SYMBOL  ( A+4. 0, G. 0, 0. 1,3, 0.0,-1 ) 

CALL  SYMBOL  (A+5. 0, G. 0, 0. 1, 3, 0.0,-1 ) 
C     THE  ABOUE  PLACES  TICK  MARKS  ALONG  THE  X  AXIS 

CALL  SYMBOL  (A+0.5, G.4, .7/G. ,  23HCUMULATIUE  POROSITY  (/),30.0,23) 

CALL  SYMBOL  ( A+0.G05, G. 0, . 7/G. 0, IHO, 0 . 0, 1 ) 

CALL  SYMBOL  ( A+0.G05, 7.0, .7/G. 0,2H50, 0.0,2) 

CALL  SYMBOL  (A+0.605, 8. 0, .7/G.0,3H100, 0.0, 3) 

CALL  SYMBOL  (A+0.G05,8.9, .7/G.O, 3H150, 0.0, 3) 
C     THE  ABOUE  PLOTS  AND  LABELS  THE  Y  AXIS 

CALL  SYMBOL  (A+1. 0, 7. 0, 0. 1,3, 0. 0, -1) 

CALL  SYMBOL  (A+1 . 0, 8. 0, 0. 1, 3, 0.0, -1 ) 

CALL  SYMBOL  (A+1 . 0, 9.0. 0. 1, 3, 0. 0,-1) 
C     THE  ABOUE  PLACES  TICK  MARKS  ALONG  THE  Y  AXIS 
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CALL  PLOT 

CALL  PLOT 

CALL  PLOT 

CALL  PLOT 

c 

THE  ABOUE 

c 

THESE  CAR; 

CALL  SYMBOL  (A+4.40. 8.80, .085,EHLEGEND. O.O.B) 
THE  ABOUE  WRITES  LEGEND  ON  THE  PLOT 
CALL  PLOT  (A+1.0.3.0.3) 

(A+G.0,9.0,2) 

(A+G. 0,6.0.2) 

(A+1.0,G,0,2) 

(A+1.0.9.0.2) 

5  CARDS  DRAW  THE  OUTLINE  OF  THE  GRAPH 
THESE  CARDS  FURTHER  LABEL  THE  GRAPH  AND  PROUIDE  THE  GRID. 
CALL  PLOT(A+1.0,8.0,3) 
CALL  PLOT(A+4.0,B.0,2) 
CALL  PLOT(A+1.0,7.0,3) 
CALL  PLOT(A+G.0,7.0,2) 
CALL  PLOT(A+5.0.G.0,3) 
CALL  PLOT(A+5.0,7.5,2) 
CALL  PLOT(A+4.0,9.0,3) 
CALL  PLOT(A+4.0,G.0,2) 
CALL  PLOT(A+3.0,G.0.3) 
CALL  PLOT(A+3.0,S.0,2) 
CALL  PLOT(A+2.0,S.0,3) 
CALL  PLOT(A+2.0,G.0,2) 
DO  2   I0U=1,5 
CI0U=I0U 
DO  1   10=2,9 
CIOA=IO 

CIO=ALOG10(CIOA) 
CALL  SYMBOL  (A+CIOU+CIO,G.O. .05,3,0.0,-1) 

1  CONTINUE 

2  CONTINUE 

3  CONTINUE 

CALL  PLOT  (A+1.0,G.0,3) 
DO  4   I=2,M+N+1 
D(I-1)=D(I) 

4  CUSP(I-1)=CUSP(I)/P0R0S 
D(M+N+l)=-3.0-A 
D(M+N+2)=1.0 
CUSP(M+N+l)=-3.00 
CUSP(M+N+2)=0.5 

CALL  SYMBOL  ( A+4. 10, 8.59-B, .07,NL, 0.0,-1) 

CALL  SYMBOL(A+4.30,8.55-B, .0S0,11HSAMPLE  CODE, 0.0, 11) 

CALL  SYMBOL(A+5.20,8.55-B,0.08,TITLE,0.0,30) 
C 

C     THE  ABOUE  LISTS  THE  WATER  CONTENT  AND  PERMEABILITY  FOR  THE  PORE 
C     SIZE  CURUE 
C 

CALL  LINE  (D. CUSP, N+N, 1,1, NL) 
C 

C   ■   THE  ABOUE  PLOTS  THE  ACTUAL  PORE  SIZE  CURUE 
C     STEPS  BELOW  PLOT  TOTAL  POROSITY  ON  THE  ORDINATE 
C 
C     B  IS  FOR  THE  PLACEMENT  OF  THE  PERMEABILITY  AND  WATER  CONTENT 

C     UALUES. 

DO  5   1=1, N+M 
IA=N+M+1-I 

5  D(IA+1)=D(IA) 
RETURN 

END 
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PROGRAM  PORE  ( INPUT, OUTPUT. TPlPE5=INPUT»TAPEE=0UTPUT, PLOT) 
C 

C    THIS  PROGRAM  PLOTS  SUCTION  US  WATER  CONTENT  GRAPHS  OBTAINED 
C    FROM  CONUENTIONAL  METHODS  AND  MERCURY  INTRUSION  TECHNIQUE 
C    IN  THE  SAME  AXES  FOR  THE  PURPOSE  OF  COMPARISON 
C 

INTEGER   TITLEO) 

REAL  DWC(30),SUCTN(30) 

REAL  PS(50),S(50).P(50),CU(50),PU1(50).PW(50).NC(50) 

A=0.0 

C=0.0 

READ.NSET 

CALL  PLOTS 

CALL  PL0T(2. 0,1.2,-3) 

DO  13  L=1,NSET 

NL=0 

E=0.0 

READ,NDATA 

READ(5,*)  W 
C 

C     W  =  SATURATED  WATER  CONTENT 
C 

PW1(1)=0. 

WC(1)=W 

c 

READ (5.*)  NUMSUC 
C     NUMSUC  =  NUMBER  OF  EXPERIMENTAL  READINGS 
C 

WRITECG.l) 

1  F0RMAT(lHl,5(/),TE5,?iC0NUENTI0NAL   METHOD?:) 
WRITE(G,2) 

2  FORMAT (/////T20,?:WATER   CONTENT?',  T40,;iSUCTION?:,/,T22,?!(PERCENT)?i, 
1  T39,^(CM   WATER):!^//) 

C 

C     READ  EXPERIMENTAL  SUCTION  AND  WATER  CONTENT  DATA 

C 

DO  4   1=1, NUMSUC 

READ(5,»)   DWC(I),SUCTN(I) 

WRITE(G,3)  DWC(I).SUCTN(I) 

3  FORMAT(T25,F5.2,T40,F8.1) 

SUCTN  ( I )  =1 .  G*AL0G1 0  ( SUCTN  CD) 

4  CONTINUE 
C 

C     DWC  =  WATER  CONTENT 
C     SUCTN  =  SUCTION 
C 

DO  12  J=1,NDATA 

READ(5,5)  TITLE 

5  FORMftT(3A10) 
C 

C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  BE  ABLE  TO  CALCULATE  RESULTS  FOR 

C     MORE  THAN  ONE  SET  OF  DATA  IN  A  SINGLE  RUN.   NDATA  IS  THE  NUMBER 

C     OF  SETS  OF  DATA  TO  BE  CALCULATED.   NDATA  IS  AN  INTEGER  UALUE,  AND 

C     IS  THE  FIRST  DATA  CARD  FOR  THE  SET  OF  DATA. 

C     DETERMINATION  OF  PORE  SIZE  DISTRIBUTION 

C 

C 

C 
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REPlD.  WS,  DS.  DM,  UP.  WSP.  WSPM.  PE,  PF,  SR 
C 

C     WS=UIEIGHT  OF  SAMPLE,  DS=SPECIFIC  GRAUITY  OF  SOLIDS,  DM=DENSITY  OF 
C     MERCURY,  UP=UOLUME  OF  PENETROMETER,  WSP=UEIGHT  OF  SAMPLE  AMD  PEM- 
C     ETROMETER,  USPM=WEIGHT  OF  SAMPLE  ,  PENETROMETER  AND  MERCURY. 
C     PE  =EUACUATING  PRESSURE,  PF=FILLING  PRESSURE, IN  MM  OF  MERCURY, 
C     SR  =STEM  READING  AT  FILLING  PRESSURE.   THE  ABOUE  UALUES  ARE  ALL  ON 
C     THE  SAME  DATA  CARD 
C 

WM=USPM  -  USP 
UM=WN/DM 
PU=PE*UM 
UAI=PU/PF 
C 

READ.N,M.SI 
C 

C     N=INTEGER  UALUE  FOR  THE  NUMBER  OF  LOW  PRESSURE  READINGS  MADE,  NOT 
C     COUNTING  THE  FILLING  PRESSURE  READING.   M=INTEGER  UALUE  FOR  THE 
C  '   NUMBER  OF  HIGH  PRESSURE  READINGS  MADE,  NOT  COUNTING  THE  INITIAL 
C     POROSINETER  READING.   SI=INTRUSION  READING  AT  INITIATION  PRESSURE 
C     IN  ML.   THE  THREE  UALUES  ARE  ALL  ON  ONE  DATA  CARD. 
C 

C     LOW  PRESSURE  CALCULATIONS 
C 

CU(l)  =  0. 
DO  8   1=2, N+1 
READ.PS(I),S(I) 
C 

C     P(I)  AND  S(I)  ARE  THE  PRESSURE  AND  INTRUSION  UALUES  FOR  EACH 
C     READING,  P(I)  IS  GIUEN  IN  MM  OF  MERCURY  FOR  THE  LOW  PRESSURE 
C     INTRUSION  AND  PSI  FOR  THE  HIGH  PRESSURE.   S(I)  IS  IN  ML.   THE  TWO 
C     UALUES  ARE  TYPED  ON  ONE  DATA  CARD,  WITH  A  CARD  FOR  EACH  READING. 
C 

P(I)=PS(I)*0. 01934 
UA=PU/PS(I) 
CUA=UAI-UA 
CSR=S(I)-SR 
SN=CSR-CUA 

IF  (SN.LT.O.)  GO  TO  G 
SG  =SN/WS 
GO  TO  7 
E     SG=0. 

SN  =0. 
7     CU(I)  =  CU(I-n  +  SG 
WC(I)=W-CU(I)*100. 
PW(I)=12.12S811*P(I) 
C 

C     12.129811  IS  THE  CONUERSION  FACTOR 
C 

PW1(I)=  l.G«ALOG10(PWCI)) 
C 

C     P(I)=INTRUSION  PRESSURE  IN  MM  OF  MERCURY,  P3=INTRUSI0N  PRESSURE  IN 
C     PSI,  UA=UOLUME  OF  AIR  IN  ML,  S(n=STEM  READING,  SN=CORRECTED 
C     INTRUSION  FOR  EACH  INCREMENT  IN  ML,  SG=CORRETED  INTRUSION  /GRAM 
C     FOR  EACH  INCREMENT,  CU( I )=CUMULATIUE  CORRECTED  INTRUSION/GRAM, 
C     D(I)=LINITING  PORE  SIZE  DIAMETER 
C 

UAI=UA 
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8  SR=S(I) 

C 

C  HIGH  PRESSURE  CftLCULPTIONS 

C 

HCI=11.».0011»UM/15000. 

DO  9   I=M+2.M+N+1 

READ.P(I),S(I) 

P(I)=P(I)+11. 
C 

C     HC  IS  CORRECT I OM  FOR  COMPRESSION  OF  MERCURY 
C 

HC=P(I)».0011»UM/15000. 

CHC=HC-HCI 

CSR=SCn-SI 

SN=CSR-CHC 

SGH=SN/-WS 

CU(I)=CU(I-1)+SGH 

WC(I)=W-CU(I)*100. 

PW(I)=P(I)*12. 123811 

PU1(I)=  l.G*ALDG10(PW(in 
C 

C     P(I)=IHTRUSION  PRESSURE  IN  PSI,  S(I)=STEM  READING,  HC=MERCURY 
C     CORRECTION,  CHC=CHANGE  IN  MERCURY  CORRECTION  IN  ML,  SN=C0RRECTEI3 
C     INTRUSION  IN  ML,  SGH=CORRECTED  INTRUSION  /GRAM,  CU( I)=CUMULATIUE 
C     CORRECTED  INTRUSION/GRAM 
C 

HCI=HC 
3     SI=S(I) 
C 

WRITE(G,10) 

10  F0RMAT(lHl,S(/),T25,;iMERCURY  INTRUSION  METHOD?!) 
WRITE(G,2) 

DO  11   I  =  2,N+M+1 

WRITE(G,3)  UCCn.PWCI) 

11  CONTINUE 
C 

CALL  JPLOTC PUl , WC, M, N, A, NL, C, J, E, TITLE, DWC, SUCTN, NUMSUC) 

NL=NL+1 

E=E+.40 

12  CONTINUE 
A=A+10.0 
C=C+50.0 

13  CONTINUE 

CALL  PLOTCO. 0,0.0,399) 
STOP 
END 
C 

C      »»»»»»«»»»»«»»»»»»»»«»»**»»-JHHt**»»*»»»»»»»1Ht«»»»»JHHt')Ht«««-)Ht«*»* 
C 

SUBROUTINE  JPLOTC PWl , UIC, M, N, A, NL, C, J, E, TITLE, DWC, SUCTN. NUMSUC ) 

INTEGER  TITLE(3) 

REAL  MC,UC(50),PU(50),Pm(50),DWC(30),SUCTN(30) 

CALL  FACTOR (0.8) 

IFCJ.GT.l)  GO  TO  21 

CALL  SYMBOL(A+3.17,0.9,.7/G.,17HWATER  CONTENT  (  ),0.0,17) 

CALL  SYMBOL (A+1., 1.25  , .7/G  ., 

1  GIHO         5        10        15        20        25        30 

2  ,0.,G1) 
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DO    14      1=2,7 
CI=I 
CALL   SYMBOL (A+CI,l.e    .    .1,13»0.    .-1    ) 

14  CONTINUE 

CALL  SYMBOL(A+  .5,4  .  , .7/6. .SSHSUCTION   (CM  OF  WATER), 90. 0,22) 

CALL  SYMBOL (A+. 7, 1.G,.7/B..1H1,0.,1) 

CALL  SYMBOL  (A+.G, 2»1.E, .7/G. ,2H1G. 0. ,2) 

CALL  SYMB0L(A+.B,3»1.G, .7/G. ,3H100, 0. . 3) 

CALL  SYMBOL(A+.G,4.0»1.G, .7/G. ,4H1000, 0.0, 4) 

CALL  SYMBOL(A+.5,5.0«1.G, .7/G. .SHIOOOO, 0. .  5) 

CALL  SYMBOL(A+.4,B.O*1.B, .7/G. ,GH100000, 0. 0,G) 

DO  15   1=1,7 

CI  =  I 

BI=1.G»CI 

CALL  SYMB0L(A+1.,BI, 0.1, 15,0.  ,-1) 

15  CONTIHUE 

CALL  PL0T(A+1.,1.G,3) 

CALL  PLOT(A+l.,G.0*l.G,2) 

CALL  PLOT(A+7.,G.0*l.G,2) 

CALL  PL0T(A+7.,1.G,2) 

CALL  PL0T(A+1. ,1.6.2) 

DO  IG   1=1,4 

Y1=FL0AT(I) 

CALL  PLOT(A+3.0+Y1, 1.6,3) 

CALL  PLOTCA+3.0+Y1,G.O»1.G,2) 

16  CONTIMUE 

CALL  PLOTCA+P.0,3.0*1.6,3) 

CALL  PL0T(A+2. 0,6.0*1.6. 2) 

CALL  PLOT(A+3.0,G.O*1.6,3) 

CALL  PLOT(A+3.0,3.0»1.6,2) 

CALL  PLOT(A+4.0,2.0«1.G,3) 

CALL  PLOT(A+7.0,2.0*1.6,2) 

CALL  PLOT(A+7.0,3.0*1.G,3) 

CALL  PLOT(A+1.0,3.0-1.6,2) 

CALL  PLOT(A+1.0,4.0»1.6,3) 

CALL  PLOT(A+7.0,4.0*1.6,2) 

CALL  PLOT(A+7.0,5.0*1.6,3) 

CALL  PLOT(A+1.0,5.0*1.6,2) 

DO  18   I0U=1,5 

CI0U=I0U 

BI0U=1.6*CI0U 

DO  17   10=2.9 

CIOA=IO 

CIO=1.G*ALOG10(CIOA) 

CALL  SYMBOL (A+1.,BIOU+CIO,0.05» 15, 0.»-l) 

17  CONTINUE 

18  CONTINUE 

CALL  PL0T(A+1., 1.6,3) 

DO  20   10=1,6 

CI0=IO 

DO  19   I0U=1,5 

CI0U=I0U 

CIGU=CI0U/5. 

CALL  SYMBOL(A+CIO+CrOU,l.G  ,.05,13,0.  ,-1) 

19  CONTINUE 

20  CONTINUE 

21  CONTINUE 

CALL  SYMBOL (A+1 . 4, 4. 04, 0 . 07, 5,0.0,-1) 
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CALL  SYMBOL(A+1.E»4.0,O.OS.19HEXPERIMENTAL  METHOD, 0.0. 19) 
CALL  SYMBOL(A+1.4,3.G4-E,0.07,NL,0.0,-1) 
CALL  SYMBOL(A+1.B,3.G-E,0.08,11HSAMPLE  CODE, 0.0, 11) 
CALL  SYMBOLC A+2. 5, 3. G-E, 0 . 08, TITLE. 0.0,30) 

WC(M+M-5)=-5.0-C 

WC(h+N-4)=5.0 

PWl(M+N-5)=-l.G 

PUI1(M+H-4)  =  1.0 


DUIC(NUMSUC+l)=-5.0-C 

DWC(hUhSUC+2)=5.0 

SUCTN(NUMSUC+1)=-1.G 

SUCTN(MUNSUC+2)=1.0 
CALL   LirSE(WC,PWl,M+N-G,l,-l,ML) 
CALL   LIMECDUCSUCTH.MUMSUC,  1,-1,5) 
RETURN 
END 
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APPENDIX  D 

Purdue  Universicy  Negative  Numbers  for  Photographs 

Figure  Negative  Number 

3  75529-24 

4  75529-34 

5  75529-25 

6  75529-28 

7  75609-14 

8  82436-11 

9  82436-5 
10  82436-2 


